BERGMAN METRICS AND GEODESICS IN THE SPACE OF KAHLER 

METRICS ON TORIC VARIETIES 



Abstract. A guiding principle in Kahler geometry is that the infinite dimensional sym- 
metric space 7i of Kahler metrics in a fixed Kahler class on a polarized projective Kahler 
manifold M should be well approximated by finite dimensional submanifolds C 7i of 
Bergman metrics of height k (Yau, Tian, Donaldson). The Bergman metric spaces are sym- 
metric spaces of type Gc/G where G = U {dk + 1) for certain dk- This article establishes the 
basic estimates for Bergman approximations for geometric families of toric Kahler manifolds. 

The approximation results are applied to the endpoint problem for geodesies of which 
are solutions of a homogeneous complex Monge- Ampere equation in A x X, where A C C 
is an annulus. Donaldson, Arezzo-Tian and Phong-Sturm raised the question whether 7i- 
geodesics with fixed endpoints can be approximated by geodesies oiBk- Phong-Sturm proved 
weak C^-convergence of Bergman to Monge- Ampere geodesies on a general Kahler manifold. 
Our approximation results show that one has C'^{A x X) convergence in the case of toric 
Kahler metrics, extending our earlier result on CP^. 

In subsequent papers, the techniques of this article are applied to approximations for 
harmonic maps into Ti., to test configuration geodesic rays and to the smooth initial value 
problem. 
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1. Introduction 

This is the first in a series of articles on the Riemannian geometry of the space 

n = {if eC°°{M) -.uj^ = ujo + dd'if > } (1) 

of Kahler metrics in the class [uq] of a polarized projective Kahler manifold {M,uq, L), 
equipped with the Riemannian metric g-j-i of Mabuchi-Semmes-Donaldson [Mlt [Si ID2j . 

r to"'' 

= / 1^1'^' where ^ en and i:eT^nc^C°^{M). (2) 

Here, L — > M is an ample line bundle with Ci(L) = [ujq]. Formally, {Ti-jgy) is an infinite di- 
mensional non-positively curved symmetric space of the type Gq/G, where G = SDif fi^g{M) 
is the group of Hamiltonian symplectic diffeomorphisms of (M, ci^o). This statement is only 
formal since G does not possess a complexification and Ti is an incomplete, infinite dimen- 
sional space. An attractive approach to the infinite dimensional geometry is to approximate 
it by a sequence of finite dimensional submanifolds Bk C 7i of so-called Bergman (or Fubini- 
Study) metrics. The space Bk of Bergman metrics may be identified with the finite dimen- 
sional symmetric space GL{dk + 1, C)/U{dk + 1) where d^ is a certain dimension. Thus, Bk is 
equipped with a finite dimensional symmetric space metric ^fg^. , which is not the same as the 
sub manifold Riemannian metric induced on it by gn- The purpose of the series is to show 
that much of the symmetric space geometry of {Bk,gBk) tends to the infinite dimensional 
symmetric space geometry of {H.gn) as /c ^ oo. 
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To put the problem and results in perspective, we recall that at the level of individual 
metrics G 7Y, there exists a well-developed approximation theory: Given a;, one can define 
a canonical sequence of Bergman metrics uJk G which approximates uj in the C°° topology 
(see (j9])), in much the same way that smooth functions can be approximated by Bernstein 
polynomials (Yau ^ and Tian [T]; see also [Cl IZelt IZe2j ). The approximation theory is 
based on microlocal analysis in the complex domain, specifically Bergman kernel asymptotics 
on and off the diagonal |BSj[ O IZell IDll IPS3] . As will be shown in jRZ3j . one may use 
the same methods to prove that the geometry of {Bk,gBk) tends to the geometry of {H.gn) 
at the infinitesimal level: e.g. that the Riemann metric, connection and curvature tensor 
of Bk tend to the Riemann metric, connection and curvature of Ti. But our principal aim 
in this series is to extend the approximation from pointwise or infinitesimal objects to more 
global aspects of the geometry, such as such as Bk -geodesies or harmonic maps to {Bk, QBk)- 
These more global approximation problems are much more difficult than the infinitesimal 
ones. The obstacles are analogous to those involved in complexifying SDif f^^lM). We will 
explain this comparison in more detail in §1.61 at the end of the introduction. 

This article is concerned with the approximation of gf-^-geodesic segments Ut in 7i with 
fixed endpoints by (y's^.-geodesic segments in Bk- As recalled in §1.11 the geodesic equation 
for the Kahler potentials (pt of Ut is a complex homogeneous Monge- Ampere equation. Little 
is known about the solutions of the Dirichlet problem at present beyond the regularity result 
that (fit £ C^'"([0,T] X M) for all a < 1 if the endpoint metrics are smooth (see X. Chen 
|Ch] and Chen-Tian |CT] for results and background). It is therefore natural to study the 
approximation of Monge- Ampere (^-n-geodesics (ft by the much simpler ^fg^. -geodesies fkit, z), 
which are defined by one parameter subgroups of GL{dk + 1, C) (see (l24l) ). The problem of 
approximating 7Y-geodesic segments between two smooth endpoints by S^-geodesic segments 
was raised by Donaldson |Dlj . Arezzo-Tian [ATj and Phong-Sturm [PSlj and was studied 
in depth by Phong-Sturm in |PSll IPS2] . They proved in |PSlj that (fkit, z) — > y9f in a weak 
C° sense on [0, 1] x M (see (IT^ ): a C° result with a remainder estimate was later proved by 
Berndtsson [B] for a somewhat different approximation. 

In this article, we study the (7^^. -approximation of (77.^-geodesics in the case of a polarized 
projective toric Kahler manifold. Our main result is that a g-H geodesic segment of toric 
Kahler metrics with fixed endpoints is approximated in by a sequence ipk{t,z) of toric 
qbi^- geodesic segments. More precisely, for any T G IR+, (pk(t,z) — > ipt{z) in C^([0,T] x M), 
generalizing the results of |SoZlj in the case of CP^. It is natural to study convergence of two 
(space-time) derivatives since the Kahler metric = ujo + dd'^ip involves two derivatives. In 
the course of the proof, we introduce methods which have many other applications to global 
approximation problems on toric Kahler manifolds, and which should also have applications 
to non-toric Kahler manifolds. 

Here, as in |SoZ2l IRZll IRZ2] . we restrict to the toric setting because, at this stage, it 
is possible to obtain much stronger results than for general Kahler manifolds and because 
it is one of the few settings where we can see clearly what is involved in the classical limit 
as ^ 00. The simplifying feature of toric Kahler manifolds is that they are completely 
integrable on both the classical and quantum level. In Riemannian terms, the submanifolds 
of toric metrics of Ti and Bk form totally geodesic fiats. Hence in the toric case, the geodesic 
equation along the fiat is linearized by the Legendre transform, with the consequence that 
there exists an explicit formula for the Monge-Ampere geodesic (ft between two smooth 
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endpoint metrics. In particular, the explicit formula shows that geodesies between smooth 
endpoints are smooth. We use this explicit solution throughout the article, starting from (129|) . 
Thus, in the toric case we only need to prove C^-convergence of the Bergman approximation. 
An analogous result on a general Kahler manifold would require an improvement on the 
known regularity results on Monge- Ampere geodesies in addition to a convergence result. 
We refer to [CTj for the state of the art on the regularity theory. 

1.1. Background. To state our results, we need some notation and background. Let L 
M"^ be an ample holomorphic line bundle over a compact complex manifold of dimension m. 
Let ujq G if^^'^^(M, Z) denote an integral Kahler form. Fixing a reference hermitian metric 
/lo on L, we may write other hermitian metrics on L as /i^ = e~'^/iO) cind then the space 
of hermitian metrics h on L with curvature (1, l)-forms Uh in the class of ujq may (by the 
dd lemma) be identified with the space Ti of relative Kahler potentials ([1]). We may then 
identify the tangent space T^H at (p E H with C°°{M). Following |Mll [Si IDlj . we define 
the Riemannian metric ([2]) on 7i. With this Riemannian metric, Ti is formally an infinite 
dimensional non-positively curved symmetric space. 

The space Bk of Bergman (or Fubini-Study) metrics of height k is defined as follows: 
Let H^{M,L^) denote the space of holomorphic sections of the kth. power L M oi L 
and let 4 + 1 = dim if°(M, L^). We let BH^{M,L^) denote the manifold of all bases 
s = {sq, . . . , Srf^} of H^{M^ L^). Given a basis, we define the Kodaira embedding 

L^:M^Cf>^\ z^[so{z),...,Sd,{z)]. (3) 

We then define a Bergman metric (or equivalently, Fubini-Study) metric of height /c to be a 
metric of the form 

K ■■= &FsY'' = -T7^, (4) 

where hps the Fubini-Study Hermitian metric on 0{1) CP'^'=. We then define 

Bu = {hs_. seBH\M,L^)}. (5) 

We use the same notation for the associated space of potentials such that hs_ = e~'^/io and 
for the associated Kahler metrics u^. We observe that with a choice of basis of H^[M, V') 
we may identify Bk with the symmetric space GL{dk + 1, C)/U{dk + 1) since GL{dk + 1, C) 
acts transitively on the set of bases, while L*hFs is unchanged if we replace the basis s by a 
unitary change of basis. 

Several further identifications are important. The first is that Bk may be identified with 
the space Ik of Hermitian inner products on H^{M,L''), the correspondence being that a 
basis is identified with an inner product for which the basis is Hermitian orthonormal. As 
in [Dll ID4] , we define maps 

Hilbk -.H^Ik, 

by the rule that a Hermitian metric h & Ti induces the inner products on H^{M, L^), 

WsWmib^h) =R \s{z)\lkdVh, (6) 

JM 
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where dVh = '\, and where R = fAt],, ^ ■ Also, h'^ denotes the induced metric on L'^. 

" m! ' Vol{M,dVh) ' 

Further, we define the identifications 

FSk -.Tk^Bu 

as follows: an inner product G = { , ) on H^{M, L'^) determines a G-orthonormal basis 
§. = of H^[M,L^) and an associated Kodaira embedding ([3]) and Bergman metric (jl]). 
Thus, 

FSk{G) = K^. (7) 

The right side is independent of the choice of /iq and the choice of orthonormal basis. As 
observed in [Dll IPSl] , FSk{G) is characterized by the fact that for any G-orthonormal basis 
{sj} of H\M,L^), we have 

= (VzGM). (8) 

i=o 

Metrics in Bk are defined by an algebro-geometric construction. By analogy with the 
approximation of real numbers by rational numbers, we say that h E T-i (or its curvature 
form uoh) has height if /i G Bk- A basic fact is that the union 

oo 

B=[jBk 

k=l 

of Bergman metrics is dense in the C°°-topology in the space Ti (see [Xl IZelj ). Indeed, 

^^^^^L^|^ = l + 0(fc-), (9) 

where the remainder is estimated in G^{M) for any r > 0; left side moreover has a complete 
asymptotic expansion (see |D3t IPSl] for precise statements). 

Now that we have defined the spaces Ti. and Bk, we can compare Monge- Ampere geodesies 
and Bergman geodesies. Geodesies of Ti satisfy the Euler-Lagrange equations for the energy 
functional determined by ([2]); see (pSj) . By [Mil B| [D2] . the geodesies of H in this metric are 
the paths hf = c'^^^Hq which satisfy the equation 

^-^|Vvi|^, = 0, (10) 

which may be interpreted as a homogeneous complex Monge-Ampere equation on A x M 
where A is an annulus [SI ID2j . 

Geodesies in B^ with respect to the symmetric space metric are given by orbits of certain 
one-parameter subgroups a\ = e*^*^ of GL{dk + 1,C). In the identification of Bk with the 
symmetric space Xk — GL{dk + 1, C)/U{dk + 1) of inner products, the 1 PS (one parameter 
subgroup) 6*^*= G GL{dk + 1) changes an orthonormal basis S*-"-* for the initial inner product 
Go to an orthonormal basis 6*"^*= • s^^-* for Gt where Gt is a geodesic of X^. Geometrically, a 
Bergman geodesic may be visualized as the path of metrics on M obtained by holomorphically 
embedding M using a basis of H^{M, L^) and then moving the embedding under the 1 PS 
subgroup e*"^* of motions of CP'^'=. The difficulty is to interpret this simple extrinsic motion 
in intrinsic terms on M. 
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In this article, we only study the endpoint problem for the geodesic equation. We assume 
given ho, hi E TL and let h{t) denote the Monge-Ampere geodesic between them. We then 
consider the geodesic Gk{t) of 1^ between G/fc(0) = Hilbk{ho) and Gfc(l) = Hilhk{hi) or 
equivalently between FSk o Hilbkiho) and FSk o Hilhkihi). Without loss of generality, 
we may assume that the change of orthonormal basis (or change of inner product) matrix 
Cfc = e^* between Hilbk{ho), Hilbk{hi) is diagonal with entries e'^'', 6'^'*'= for some \j E M. 
Let s^*-* = e*^*-' ■ where e*^*" is diagonal with entries e^^*. Define 

hk{t) := FSk o Gk{t) = =: hoe-^'^'\ (11) 

It follows immediately from that 

We emphasize that ipk{t; z) is the intrinsic Bu geodesic between the endpoints FSkoHilhk{ho) 
and FSkoHilbk{hi). It is of course quite distinct from the Hilbk-iTaage of the Monge-Ampere 
geodesic; the latter is not intrinsic to Bk and one cannot gain any information on the Ti- 
geodesic by studying it. 

Let us summarize the notation for hermitian metrics and geodesies of metrics: 

• For any metric h on L, h^ denotes the induced metric on L^, and for any metric H 
on L'^, i/fc is the induced metric on L; 

• Given h^ eH, ht = e~'^*/io is the Monge-Ampere geodesic; 

• hk = FS o Hilbk{h) E Bk is the natural approximating Bergman metric to h, and 
hk{t) = e~'^'=^*-'/io is the Bergman geodesic ( ITT]) . 

The main result of Phong-Sturm [PSlj is that the Monge-Ampere geodesic (ft is approxi- 
mated by the IPS Bergman geodesic (pk{t, z) in the following weak C° sense: 



ipt{z) = lim 



snp(pk{t,z) 



k>e 



uniformly as £ ^ cxd, (13) 



where u* is the upper envelope of u, i.e., u*{(q) = linie^o sup|^_^jj|<^ m(C). In particular, 
without taking the upper envelope, supk>e'^k{t, z) —>■ (p{t,z) almost everywhere as i —>■ oo. 
See also [B] for the subsequent proof of an analogous result for the adjoint bundle L'' ® K 
(where K is the canonical bundle) with an error estimate ||v5fc(^) ~ V^(^)||co = 0{)^^). 

1.2. Statement of results. Our purpose is to show that the degree of convergence of 
hk{t) ht or equivalently of (pk{ii z) ^t{z) is much stronger for toric hermitian metrics 
on the invariant line bundle L ^ M over a smooth toric Kahler manifold. We recall that a 
toric variety M of dimension m carries the holomorphic action of a complex torus (C*)™ with 
an open dense orbit. The associated real torus T™ = (S*^)™ acts on M in a Hamiltonian 
fashion with respect to any invariant Kahler metric u;, i.e., it possesses a moment map 
fi : M P with image a convex lattice polytope. Here, and henceforth, P denotes the 
closed polytope; its interior is denoted P° (see §2] for background). Objects associated to M 
are called toric if they are invariant or equivariant with respect to the torus action (real or 
complex, depending on the context). We define the space of toric Hermitian metrics by 

Ht- = {y^En: (e*^)V = V, for all e'^ E T™}. (14) 
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Here, we assume the reference metric ho is T'"-invariant. We note that since T"^ has a 
moment map, it automatically lifts to L and hence it makes sense to say that ho : L C 
is invariant under it. With a slight abuse of notation carried over from [Dl] . we also let (f 
denote the full Kahler potential on the open orbit, i.e., = dd'^ip on the open orbit. It is 
clearly T™-invariant. 
Our main result is 

Theorem 1.1. Let L ^ M be a very ample toric line bundle over a smooth compact toric 
variety M. Let Tix denote the space of toric Hermitian metrics on L. Let ho, hi G Tir 
and let ht be the Monge-Ampere geodesic between them. Let hk{t) be the Bergman geodesic 
between Hilhk{ho) and Hilbk{hi) in Bk- Let hk{t) = e~'^'=(*'^)/zo and let ht = e~'^^^^^ho. Then 

lim ipk{t,z) = Lpt{z) 

in C^([0, 1] X M). In fact, there exists C independent of k such that 

\Wk- ^\\c-^{[o,i]>,M) <Ck-^'^+\ Ve>0. 

Our methods show moreover that away from the divisor at infinity T) (cf. the function 
(pk{t,z) has an asymptotic expansion in powers of k~^, and converges in C°° to (ft- But 
the asymptotics become complicated near V, and require a 'multi-scale' analysis involving 
distance to boundary facets. It is therefore not clear whether (fk has an asymptotic expansion 
in k~^ globally on M. At least, no such asymptotics follow from the known Bergman kernel 
asymptotics, on or off the diagonal. The analysis of these regimes for general toric varieties 
seems to be fundamental in 'quantum mechanical approximations' on toric varieties. 

As mentioned above, the Monge-Ampere equation can be linearized in the toric case and 
solved explicitly (IT7I) : we give a simple new proof in §2J The geodesic arcs are easily seen to 
be C°° when the endpoints are C°°. Hence the C^-convergence result does not improve the 
known regularity results on Monge-Ampere geodesies of toric metrics, but pertain only to the 
degree of convergence of Bergman to Monge-Ampere geodesies in a setting where the latter 
are known to be smooth; it is possible that the methods can be developed to give regularity 
results, but this is a distant prospect (see the remarks at the end of this introduction). 

1.3. Outline of the proof. Let us now outline the proof of Theorem II. 1[ We start 
with the fact that the Legendre transform of the Kahler potential linearizes the Monge- 
Ampere equation (cf. §2.71 and [A | O [D3] ). The Legendre transform dp of the open-orbit 
Kahler potential (f, a convex function on R'" in logarithmic coordinates, is the so-called dual 
symplectic potential 

u^{x) = C(f{x), (15) 
a convex function on the convex polytope P. Under this Legendre transform, the complex 
Monge-Ampere equation on TYt™ linearizes to the equation il = and is thus solved by 

Ut = U^o + ti^^Pi - ^vo)- (16) 

Hence the solution ipt of the geodesic equation on Ti is solved in the toric setting by 

= (17) 

Our goal is to show that ipk{t] z) C~^ut as in f|T6|) in a strong sense. 

The second simplifying feature of the toric setting occurs on the quantum level. The 
Bergman geodesic is obtained by applying the FSk map to the one-parameter subgroup 
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g*^fc general, it is difficult to understand what kind of asymptotic behavior is possessed 
by the operators e*^*. But on a toric variety, there exists a natural basis of the space of 
holomorphic sections H^{M^L^) furnished by monomial sections 2" which are orthogonal 
with respect to all torus-invariant inner products, and with respect to which all change of 
basis operators e^^'' are diagonal; we refer to ^ or to |STZ1] for background. Hence, we 
only need to analyze the eigenvalues of e^*. The exponents a of the monomials are lattice 
points a G kP in the /cth dilate of the polytope P corresponding to M. The eigenvalues in 
the toric case are given by 

1 , / 

K ■= 77 loj 



2 "\^Q,.(a 

where Q/ife(a) is a 'norming constant' for a toric inner product. By a norming constant for 
a toric Hermitian inner product G on H^[M, L^) we mean the associated norm-squares 
of the monomials 

QG{a) = \K\\l. (19) 
In particular, if /i G TYt™, the norming constants for Hilbk{h) are given by 



Qh^ia) = \K\\U ■■= / K{z)\i,dVh. (20) 



M 



p 



Thus, an orthonormal basis of H^{M, L'') with respect to Hilhk{h) for h G Tix is given by 
{ , ^" a G kP n Z™}. An equivalent, and in a sense dual (cf. [|3]), formulation is in 

terms of the functions 

12 



n.(«,-):=^^, (21) 



and their special values 



V,.{a) :=n.(a,/., (-)) = ^^^^^^ . (22) 

Given two toric hermitian metrics h^^hi G TYt™? the change of basis matrix e^'' = (Jho,hi,k 
from the monomial orthonormal basis for Hilbkiho) to that for Hilhk{hi) is diagonal, and 
the eigenvalues are given by 

5p(e^'=e^^) := {e^^"^'^) = -^^4^, a G A;P}. (23) 

Qhu[a) 

Hence, for a iS^-geodesic, ( 1T2|) becomes 

(^fe(t,z) = ilogZfc(t,2) (24) 

where 
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It is interesting to observe that the relative Kahler potential (12^ is the logarithm of an 
exponential sum, hence has the form of a free energy of a statistical mechanical problem 



with states parameterized by a G kP and with Boltzmann weights 
Thus, our goal is to prove that 



-ft 







1 , ^ f 



Sr,{Z 



k 



1.4. Heuristic proof. Let us next sketch a heuristic proof which makes the pointwise con- 
vergence obvious. The first step is to obtain good asymptotics of the norming constants 
( !20|) . As in [SoZlj . they may be expressed in terms of the symplectic potential by 

QhH{a) = [ e-'^K(^-)+<f-'^"^W))cia; (27) 



Jp 

As ^ oo the integral is dominated by the unique point x = ^ where the 'phase function' 
is maximized. The Hessian is always non-degenerate and by complex stationary phase we 
obtain the asymptotics 

The complex stationary phase (or steepest descent) method does not apply near the boundary 
dP, causing serious complications, but in this heuristic sketch we ignore this aspect. 
If we then replace each term in by its asymptotics, we obtain 

ipk{t,ef^^) r^^log g2fe(«oW+t(ni(a)-«oH)+(p,a»_ ,^28) 

The exponent {uo{a) + t{ui{a) — uo{a)) + (p, a)) is convex and therefore has a unique min- 
imum point. This suggests applying a discrete analogue of complex stationary phase to the 
sum (1281) . a Dedekind-Riemann sum which is asymptotic to the integral 

^2k{uo{a)+t{ui{a)-uo{a))+{p,a))^^ 

Taking ^ log of the integral and applying complex stationary phase gives the asymptote 

max{Mo(Q;) + t{ui{a) — Uo(a)) + (p, a)}- 

But this is the Legendre transform of the ray of symplectic potentials 

Mipo(a) + t{u^^{a) - u^^{a)), 

and thus is the Monge-Ampere geodesic. 

This is the core idea of the proof. We now give the rigorous version. 
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1.5. Outline of the rigorous proof. The main difficuhy in the proof of Theorem 11.11 
is that the norms have very different asymptotic regimes according to the position of the 
normahzed lattice point ^ relative to the boundary dP of the polytope. Even in the simplest 
case of CP'", the different positions correspond to the regimes of the central limit theorem, 
large deviations theorems and Poisson law of rare events for multi-nomial coefficients. In 
determining the asymptotics of (12^ . we face the difficulty that these Boltzmann weights 
might be exponentially growing or decaying in /c as /c ^ oo. 

To simplify the comparison between the Bergman and Monge-Ampere geodesies, we take 
advantage of the explicit solution (fTTIl of geodesic equation to re- write Zk{t, z) in the form 



where as usual ht = e '^^Hq (with (ft as in (fT7|) ). and where 

^kit, a) := 1 (30) 

One of the key ideas is that 7lk{t, a) is to at least one order a semi-classical symbol in k, i.e., 
has at least to some extent an asymptotic expansion in powers of k. Once this is established, 
it is possible to prove that 

^log 'J^k{t,a)Vh^{a,z) ^0 (31) 

in the C^-topology on [0, 1] x M. 

The proof of Theorem 11.11 consists of four main ingredients: 

• The Localization Lemma [1.21 which states that the sum over a localizes to a ball of 
radius 0{k~2+^^ around the point fit{z). Here and hereafter, S can be taken to be 
any sufficiently small positive constant. 

• Bergman/Szego asymptotics (see §4.2p . which allow one to make comparisons between 
the sum in Zk and sums with known asymptotics. 

• The Regularity Lemma [L3l which states that the summands TZk{t, a) one is averaging 
have sufficiently smooth asymptotics as /c — oo, allowing one to Taylor expand to 
order at least one around the point fit{z). 

• Joint asymptotics of the Fourier coefficients and particularly their special val- 
ues Vhk{a) in the parameters k and distance to dP (see Proposition 16. ip . We use 
a complex stationary phase method in the 'interior region' far from dP and local 
Bargmann-Fock models near dP. 

The Localization Lemma is needed not just for TZk{t,a) but also for summands which 
arise from differentiation with respect to {t,z): 

Lemma 1.2. (Localization of Sums) Let Bk{t, a) : Z'" fl kP ^ C be a family of lattice point 
functions satisfying \Bk{t,a)\ < C^k^'^ for some Co,M > 0. Then, there exists C > so 
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that for any S > 0, 



The proof is an integration by parts argument. One could localize to the smaller scale 
1^ — < C*^^^ but then the argument only brings errors of the order (log/c)"*^ for all 

M and that complicates later applications. 

The regularity Lemma concerns the behavior of the 'Fourier multiplier' Rk{t, a) fl30|) . The 
sum ( !25ll formally resembles the Berezin covariant symbol of a Toeplitz Fourier multiplier, 
i.e., the restriction to the diagonal of the Schwartz kernel of the operator; we refer to ( jSTZ2l 
IZe2] ) for discussion of such Toeplitz Fourier multipliers operators on toric varieties and their 
Berezin symbols. However, the resemblance is a priori just formal - it is not obvious that 
Rk{t,a) has asymptotics in k. As mentioned above, the nature of the asymptotics is most 
difficult near dP; it is not obvious that smooth convergence holds along V, the divisor at 
infinity. 

The purpose of introducing Rk{t,a) is explained by the following result. First, we make 
the 



Definition: We define the metric volume ratio to be the function on [0, 1] x P defined by 



det V^tix) ^^^^ 



(det V'^Uo{x)y~'{det V^Ui{x)y 



Lemma 1.3. (Regularity) The volume ratio 7loo{t,x) E C°°([0, 1] x P). Further, forO < j < 
2, 

(|y7^,(^,«) = (|y7^oo(t,^) + o(fc-^), 

where the O symbol is uniform in {t,a). 

This Lemma is the subtlest part of the analysis. If the TZk function were replaced by a 
fixed function f{x) evaluated at ^ then the convergence problem reduces to generalizations 
of convergence of Bernstein polynomial approximations to smooth functions |Ze2j . and only 
requires now standard Bergman kernel asymptotics. However, the actual Rk{t,a) do not 
apriori have this form, and much more is required for their analysis than asymptotics (on 
and off diagonal) of Bergman kernels. The analysis uses a mixture of complex stationary 
phase arguments in directions where ^ is 'not too close' to dP, while directions 'close to' dP 
we use an approximation by the 'linear' Bargmann-Fock model (see §2.61 and §6.4p . 

The somewhat unexpected k~^^^ remainder estimate has its origin in this mixture of 
complex stationary phase and Bargmann-Fock asymptotics. Both methods are valid for 
k satisfying '"^"^^ < 6k < ^'^/khgk' ^^^^ region, the stationary phase remainder is of 
order {k6k)^^ while the Bargmann-Fock remainder is of order kS^; the two remainders agree 
when 5k = A;~3, and then the remainder is 0{k~^^^). For smaller Sk the Bargmann-Fock 
approximation is more accurate and for larger 6k the stationary phase approximation is 
more accurate. This matter is discussed in detail in §6.4[ 
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The rest of the proof of the C^-convergence may be roughly outhned as follows: We 
calculate two logarithmic derivatives of e~'"-'"^^^^ Zk{t, z) of (!29l) with respect to (t,p). Using 
the Localization Lemma 11.21 we can drop the terms in the resulting sums corresponding to 
a for which |^ — Ht{z)\ > k~2~^^. In the remaining terms we use the Regularity Lemma [1.31 
to approximate the summands by their Taylor expansions to order one around Ht{z)- This 
reduces the expressions to derivatives of the diagonal Szego kernel 

for the metric /i^ on H^[M, V') induced by Monge- Ampere geodesic ht. Here, we use the 
smoothness of ht. The known asymptotic expansion of this kernel ( §4.21) implies the C^- 
convergence of e'^'^*^^^ Zk{t, z). 

As indicated in this sketch, the key problem is to analyze the joint asymptotics of norming 
constants Q\{(y) and the dual constants Vhk{a) fl22|) in {k,a). Norming constants are a 
complete set of invariants of toric Kahler metrics. Initial results (but not joint asymptotics 
in the boundary regime) were obtained in jSTZl] : norms are also an important component 
of Donaldson's numerical analysis of canonical metrics [D4] on toric varieties. In |SoZlj the 
joint asymptotics of Q^(a) were studied up to the boundary of the polytope [0, 1] associated 
to CP^. In this article, we emphasize the dual constants 



1.6. Bergman approximation and complexification. Having described our methods 
and results, we return to the discussion of their relation to Kahler quantization and to the 
obstacles in complexifying Dif f^^^{M). Further discussion is given in [RZ2] . 

We may distinguish two intuitive ideas as to the nature of Monge- Ampere geodesies. The 
first heuristic idea, due to Semmes [S] and Donaldson jDl] . is to view HCMA geodesies 
as one parameter subgroups of Gc where G = SDif f^^{M). One parameter subgroups of 
SDif f^^{M) are defined by Hamiltonian flows of initial Hamiltonians 0o with respect to 
ujq. a complexified one parameter subgroup is the analytic continuation in time of such 
a Hamiltonian flow [S], IDlj . This idea is heuristic inasmuch as Hamiltonian flows need not 
possess analytic continuations in time; moreover, no genuine complexification of SDif f^^^{M) 
exists. 

The second intuitive idea (backed up by the results of |PSlj and this article) is to view 
HCMA geodesies as classical limits of Bk geodesies. The latter have a very simple extrinsic 
interpretation as one parameter motions e*"^'=is(M) of a holomorphic embedding : M — s> 
CP'^'=. But the passage to the classical limit is quite non-standard from the point of view of 
Kahler quantization. The problem is that the approximating one parameter subgroups e*^* 
of operators on H^{M, L^), which change an orthonormal basis for an initial inner product 
to a path of orthonormal bases for the geodesic of inner products, are not apriori complex 
Fourier integral operators or any known kind of quantization of classical dynamics. 

The heuristic view taken in this article and series is that e*"^* should be approximately the 
analytic continuation of the Kahler quantization of a classical Hamiltonian flow. To explain 
this, let us recall the basic ideas of Kahler quantization. 

Traditionally, Kahler quantization refers to the quantization of a polarized Kahler man- 
ifold {M,u},L) by Hilbert spaces H^{M.,L^) of holomorphic sections of high powers of a 
holomorphic line bundle L ^ M with Chern class ci(L) = [uj\. The Kahler form determines 
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a Hermitian metric h such that Ric{h) = u. The Hermitian metric induces inner products 
Hilhk{h) on H^{M, L^). In this quantization theory, functions H on M are quantized as Her- 
mitian (Toephtz) operators H := UhkHUfik on H^{M, L^), and canonical transformations 
of (M, u) are quantized as unitary operators on H^{M, L^). Quantum dynamics is given by 
unitary groups e^^'^^ (see |BBSt [BSj| IZel] for references). 



In the case of Bergman geodesies with fixed endpoints, H should be ipQ, the initial tangent 
vector to the HCMA geodesic with the fixed endpoints. The quantization of the Hamiltonian 
flow of if should then be e^^^^ and its analytic continuation should be e*^^ . The change of 
basis operator 6*^*= should then be approximately the same as e*^^. But proving this and 
taking the classical limit is necessarily non-standard when the classical analytic continuation 
of the Hamiltonian flow of (p does not exist. Moreover, we only know that G C^'^ . 

This picture of the Bergman approximation to HCMA geodesies is validated in this article 
in the case of the Dirichlet problem on projective toric Kahler manifolds. It is verified for 
the initial value problem on toric Kahler manifolds in |RZ2j . In work in progress, we are 
investigating the same principle for general Kahler metrics on Riemann surfaces [RZ4j . 

1.7. Final remarks and further results and problems. An obvious question within 
the toric setting is whether (pk{i) — > </3t in a stronger topology than on a toric variety. It 
seems possible that the methods of this paper could be extended to C'^-convergence. The 
methods of this paper easily imply convergence for all k away from dP or equivalently the 
divisor at infinity, but the degree of convergence along this set has yet to be investigated. As 
mentioned above, we do not see why ipk should have an asymptotic expansion in fc, but this 
aspect may deserve further exploration. We also mention that our methods can be extended 
to prove C^-convergence of Berndtsson's approximations in [B] . 

In subsequent articles on the toric case, we build on the methods introduced here to prove 
convergence theorems for other geodesies and for general harmonic maps |RZlj (including 
the Wess-Zumino-Witten equation). In |SoZ2j . we develop the methods of this article to 
prove that the geodesic rays constructed in |PS2j from test configurations are C^'^ and no 
better on a toric variety. Test configuration geodesic rays are solutions of a kind of initial 
value problem; we refer to the articles [PS21 ISoZ2] for the definitions and results. For test 
configuration geodesies, the analogue of TZk is not even smooth in t. The smooth initial value 
problem is studied in |RZ2] . In a different direction, one of the authors and Y. Rubinstein 
prove a convergence result for completely general harmonic maps of Riemannian manifolds 
with boundary into toric varieties (see [Rl IRZlj ). This includes the Wess-Zumino-Witten 
model where the manifold is a Riemann surface with boundary. 

We believe that the techniques of this paper extend to other settings with a high degree 
of symmetry, such as Abelian varieties and other settings discussed in |D5j . The general 
Kahler case involves significant further obstacles. A basic problem in generalizing the results 
is to construct a useful localized basis of sections on a general (M, cj). In the toric case, 
we use the basis of T™'-invariant states Sa = z"--, which 'localize' on the so-called 'Bohr- 
Sommerfeld tori', i.e. the inverse images /i~^(^) of lattice points under the moment map 
Such Bohr-Sommerfeld states also exist on any Riemann surface; in |RZ4] . we relate them 
to the convergence problem for HCMA geodesies on Riemann surfaces. 

We briefly speculate on the higher dimensional general Kahler case. There are a number of 
plausible substitutes for the Bohr-Sommerfeld basis on a general Kahler manifold. A rather 
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traditional one is to study the asymptotics of e^^ on a basis of coherent states $]^fe. Here, 
^^kiz) = '^hk'y^''^) ^ g^j^g ^2 normalized Szesro kernels pinned down in the second argument. 

Intuitively, is like a Gaussian bump centered at w with shape determined by the metric 
h. It is thus more localized than the monomials z", which are only Gaussian transverse to 
the tori. Under the change of basis operators e*"^*, both the center and shape should change. 
Like the monomials z", coherent states have some degree of orthogonality. There are in 
addition other well-localized bases depending on the Kahler metric which may be used in 
the analysis. 

Our main result (Theorem II. ip may be viewed heuristically as showing that as A; — *• oo 
the change of basis operators e^^^ tend to a path ft of diffeomorphisms changing the initial 
Kahler metric lo^ into the metric lOt along the Monge- Ampere geodesic. We conjecture that 
gt^fc<|)W'^ ~ where ht is the Monge- Ampere geodesic and ft is the Moser path of 

diffeomorphisms such that fl^Q = ujf We leave the exact degree of asymptotic similarity 
vague at this time since even the regularity of the Moser path is currently an open problem. 

Acknolwedgements The authors would like to thank D. H. Phong and J. Sturm for their 
support of this project, and to thank them and Y. A. Rubinstein for many detailed correc- 
tions. The second author's collaboration with Y. A. Rubinstein subsequent to the initial 
version of this article has led to a deepened understanding of the global approximation 
problem, which is reflected in the revised version of the introduction. 



In this section, we review the necessary background on toric Kahler manifolds. In addition 
to standard material on Kahler and symplectic potentials, moment maps and polytopes, we 
also present some rather non-standard material on almost analytic extensions of Kahler 
potentials and moment maps that are needed later on. We also give a simple proof that the 
Legendre transform from Kahler potentials to symplectic potentials linearizes the Monge- 
Ampere equation. 

Let M be a complex manifold. We use the following standard notation: ^ = — 
'^'§y)iM ~ K^"'"^^)- often find it convenient to use the real operators d = d + d^dF := 
±{d-d) and dd' = j^dd. 

Let L ^ M be a holomorphic line bundle. The Chern form of a Hermitian metric h on L 
is defined by 



where denotes a local holomorphic frame (= nonvanishing section) of L over an open 
set U C M, and WeiWh = h(eL,eLY^'^ denotes the /i-norm of cl- We say that {L,h) is 
positive if the (real) 2-form Uh is a positive (1,1) form, i.e., defines a Kahler metric. We 
write II 6^(2;) 11^ = e"*^ or locally h = e""^, and then refer to ip as the Kahler potential of Uh 
in U. In this notation. 



2. Background on toric varieties 




(33) 




(34) 
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If we fix a Hermitian metric Hq and let h = e~'^ho, and put ooq = ooho, then 

Uh = ujq + d(fip. (35) 

The metric h induces Hermitian metrics on = L ® ■ ■ ■ ® L given by Hs^'^H/jj^ = 

We now specialize to toric Kahler manifolds; for background, we refer to pilD3l[GllSTZlj . 
A toric Kahler manifold is a Kahler manifold (M, J, cu) on which the complex torus (C*)'" 
acts holomorphically with an open orbit M°. Choosing a basepoint mo on the open orbit 
identifies M° = (C*)™ and give the point z = e^/^"'"*'^mo the holomorphic coordinates 

z = 6^/2+^^ e (C*)"*, p, ^ e M™. (36) 

The real torus T"* C (C*)™ acts in a Hamiltonian fashion with respect to u. Its moment 
map yU. = /i^:M— s>Pct*~ M™ (where t is the Lie algebra of T*") with respect to u 
defines a singular torus fibration over a convex lattice polytope P; as in the introduction, P 
is understood to be the closed polytope. We recall that the moment map of a Hamiltonian 
torus action with respect to a symplectic form u is the map : M ^ t* defined by 
d{fiuj{z), = L^#uj where i^* is the vector field on M induced by the vector ,^ G t. Over the 
open orbit one thus has a symplectic identification 

We let X denote the Euclidean coordinates on P. The components (Ji, . . . , J^) of the moment 
map are called action variables for the torus action. The symplectically dual variables on T"^ 
are called the angle variables. Given a basis of t or equivalently of the action variables, we 
denote by {^} the corresponding generators (Hamiltonian vector fields) of the T"^ action. 

Under the complex structure J, we also obtain generators of the action. 

The action variables are globally defined smooth functions but fail to be coordinates at 
points where the generators of the T™ action vanish. We denote the set of such points by 
V and refer to it as the divisor at infinity, li p ^ V and denotes the isotropy group of 
p, then the generating vector fields of become linearly dependent at P. Since we are 
proving estimates, we need to replace them near points of V by vector fields with norms 
bounded below. We discuss good choices of coordinates near points of V below. 

We assume M is smooth and that P is a Delzant polytope. It is defined by a set of linear 
inequalities 

£r{x) := (x, Vr) — Aj. > 0, r = 1, ...,d, 

where Vr is a primitive element of the lattice and inward-pointing normal to the r-th (m — 1)- 
dimensional facet Fj. = {ir = 0} of P. We recall that a facet is a highest dimensional face 
of a polytope. The inverse image fi~^{dP) of the boundary of P is the divisor at infinity 
P C M. For X e dP we denote by 

J^(x) = {r : 4(x) = 0} 

the set of facets containing x. To measure when a; G P is near the boundary we further 
define 

T,{x) = {r:\ir{x)\<e}. (37) 

The simplest toric varieties are linear Kahler manifolds (V, uj) carrying a linear holomorphic 
torus action. They provide local models near a corner of P or equivalently near a fixed point 
of the T*" action. As discussed in |GSt ILT] . a linear symplectic torus action is determined 
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by a choice of m elements (3j of the weight lattice of the Lie algebra of the torus. The vector 
space then decomposes (y,uj) = 0(1^, cui) of orthogonal symplectic subspaces so that the 
moment map has the form 

flBpivi, ■ ■ ■ ,Vm) = ^\Vj\'^Pj. (38) 

The image of the moment map is the orthant M™. This provides a useful local model at 
corners. We refer to these as Bargmann-Fock models; they play a fundamental role in this 
article (cf. ^23]). 

2.1. Slice-orbit coordinates. We will also need local models at points near codimension 
r faces, and therefore supplement the coordinates (136!) on the open orbit with holomorphic 
coordinates valid in neighborhoods of points of V. An atlas of coordinate charts for M 
generalizing the usual affine charts of CP"^ is given in |STZ1] . §3.2 and we briefly recall the 
definitions. For each vertex vq G P, we define the chart f/^,, by 

U,,:={zeMp;Xvo{z)y^O}, (39) 

where 

x„(^) = ^" = zr---^:". 

Throughout the article we use standard multi- index notation, and put |a| =«! + ■■■ + a„,. 
Since P is Delzant, we can choose lattice points a^, . . . , in P such that each is in an 
edge incident to the vertex wo? and the vectors :=«■' — form a basis of . We define 

: (C*)- - (C*)-, ri{z)=ri,{z):={z^\...,z^"'). (40) 

The map r/ is a T'^-equivariant biholomorphism with inverse 

z:(CX^(C*)™, z{ri) = {rf^\...,rf^"'), (41) 

where e-' is the standard basis for C", and F is an m x m-matrix with det F = ±1 and integer 
coefficients defined by 

Ff-' = 6"', = — vq. (42) 

The corner of P at t'o is transformed to the standard corner of the orthant by the affine 
linear transformation 

f : 9 n ^ Fn - Fi;o G W\ (43) 

which preserves Z™, carries P to a polytope C {x G ; > 0} and carries the facets 
Fj incident at to the coordinate hyperplanes = {x G Q^^ ; = 0}. The map r/ extends a 
homeomorphism : 

T] : f/^,(j C™", = 0, Zq = the fixed point corresponding to v^. (44) 

By this homeomorphism, the set fip^{Fj) corresponds to the set {rj G ; rjj = 0}. If F be a 
closed face with dim F = m — r which contains vq, then there are facets Fj^, . . . , Fi^ incident 
at vq such that F = Fj^ fl ■ ■ ■ fl Fj^. The subvariety ii~p'{F) corresponding F is expressed by 

/ipi(F)nf/„„ = {r/GC'";r/,^=0, j = l,...,r}. (45) 

When working near a point of /ip^(F), we simplify notation by writing 

= [r]\ 7]") G = C" X C™-" (46) 
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where 77' = {rji.) as in ( H5|) and where rj" are the remaining ?7j's, so that (0,77") is a local 
coordinate of the submanifold fip^{F). When the point (0,77") lies in the open orbit of 
fj,p^(F), we often write rj" = e*^"+^"/2. in practice, we simplify notation by tacitly treating 
the corner at vq as if it were the standard corner of MJ^, omit mention of F and always use 
{z', z") instead of rj. It is straightforward to rewrite all the expressions we use in terms of 
the more careful coordinate charts just mentioned. 

These coordinates may be described more geometrically as slice-orbit coordinates. Let 
Pq G fXp^{F) and let (C*)^^ denote its stabilizer (isotropy) subgroup. Then there always 
exists a local slice at Pq, i.e., a local analytic subspace S <Z M such that Pq & S, S is 
invariant under (C*)pg, and such that the natural (C*)™ equivariant map of the normal 
bundle of the orbit (C*)"* ■ Pq, 

[C,P]G(C*)"^X(c.)j.^^C-^eM (47) 

is biholomorphism onto {'C*)™'-S. The terminology is taken from (see Theorem 1.23). The 
slice 5* can be taken to be the image of a ball in the hermitian normal space Tpg((C*)'"Po)"'" 
to the orbit under any local holomorphic embedding w : Tpg((C*)'"Po)"'" M with w{Pq) = 
Po,dwpg = Id. The affine coordinates rj" above define the slice S = ri~^{(z' , z"{Po)) : z' G 
(C*)*"}. The local 'orbit-slice' coordinates are then defined by 

P = (z', e^^"+^"/2) ^ r]{P) = e''"+p"/\z', 0) (48) 

where {z', 0) G S* is the point on the slice with affine holomorphic coordinates z' = {rj'). 

As will be seen below, toric functions are smooth functions of the variables e''^ away from 
V, and of the variables {zjl"^ at points near V. We introduce the following 'polar coordinates' 
centered at a point P G P: 

rr.= \zj\ = e^^/'. (49) 
They are polar coordinates along the slice. The gradient vector field of rj is denoted 

As with polar vector fields, it is not well-defined at rj = 0. But to prove estimates of 
functions which are smooth functions of r| it is sufficient to prove estimates with respect 
to the vector fields 7^ or ttt^. 

drj d{Tj) 

2.2. Kahler potential in the open orbit and symplectic potential. Now consider the 
Kahler metrics uj inTi (cf. ([I])). We recall that on any simply connected open set, a Kahler 
metric may be locally expressed as = 2idd(p where (p is a locally defined function which 
is unique up to the addition (p ip + f{z) + f{z) of the real part of a holomorphic or 
antiholomorphic function /. Here, a G M is a real constant which depends on the choice of 
coordinates. Thus, a Kahler metric uj E Ti has a Kahler potential p over the open orbit 
M° C M. In fact, there is a canonical choice of the open-orbit Kahler potential once one 
fixes the image P of the moment map: 

<^(z) =Iog^|2°|2 = log^e<°'''\ (50) 

Invariance under the real torus action implies that p only depends on the p-variables, so 
that we may write it in the form 



p{z)=p{p) = F{eP). 



(51) 
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The notation (f{z) = (f{p) is an abuse of notation, but is rather standard since [D3] . For 
instance, the Fubini-Study Kahler potential is ip{z) = log(l + = log(l + e^) = F^e^). 
Note that the Kahler potential log(l + extends to C™ from the open orbit (C*)™, 
although the coordinates (p, 9) are only valid on the open orbit. This is a typical situation. 
On the open orbit, we then have 

2 ^ dpkdpj Z.J Zk 

Positivity of uJ^p implies that = F{e^) is a strictly convex function of p G M". The 

moment map with respect to uJ^p is given on the open orbit by 

p^^(^i, . . . , = V,^(p) = V,F{eP\. . . , e"-), (z = e"/'^''). (53) 

Here, and henceforth, we subscript moments maps either by the Hermitian metric h or by 
a local Kahler potential Lp. The formula (!53l) follows from the fact that the generators ^ 

of the T"^ actions are Hamiltonian vector fields with respect to with Hamiltonians 
since 

, 9^ 



L_a_uj^ = d— — . (54) 



" dpj 



The moment map is a homeomorphism from p G M™" to the interior P° of P and extends 
as a smooth map from M ^ P with critical points on the divisor at infinity V. Hence, the 
Hamiltonians fl5^ extend to P. 

Note that the local Kahler potential on the open orbit is not the same as the global smooth 
relative Kahler potential in ([1]) with respect to a background Kahler metric uq. That is, given 
a reference metric lvq with Kahler potential ipo, it follows by the dd lemma that uj = ujQ + dd^ip 
with ip G C°°{M). As discussed in |D3j (see Proposition 3.1.7), the Kahler potential ip on 
the open orbit defines a singular potential on M which satisfies dd'^ip = u + H where H is 
a fixed current supported on V. We generally denote Kahler potentials by ip and in each 
context explain which type we mean. 

By ( !52|) . a T'^-invariant Kahler potential defines a real convex function on p G M™. Its 
Legendre dual is the symplectic potential u^: for x G P there is a unique p such that 
l^tfii^^^'^) = ^p'P = X. Then the Legendre transform is defined to be the convex function 

u^{x) = (x,p^.) - V5(px), e''"/^ = p;;^(x) p^ = 21ogp^^(x) (55) 

on P. The gradient VxU^ is an inverse to on on the open orbit, or equivalently on 
P, in the sense that Vuy,{p^^{z)) = z as long as Puj^{z) ^ dP. 

The symplectic potential has canonical logarithmic singularities on dP. According to 
(Proposition 2.8) or [D3j ( Proposition 3.1.7), there is a one-to-one correspondence between 
T^-invariant Kahler potentials ip on Mp and symplectic potentials u in the class S of 
continuous convex functions on P such that u — uq is smooth on P where 



Uo[X) 



^4(x)log4(x). (56) 



Thus, u^{x) = uo{x) + fip{x) where G C°°{P). We note that uq and are convex, that 
Mo = on dP and hence = on dP. By convexity, maxp uq = 0. 
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We denote by = V^m^ the Hessian of the symplectic potential It has simple poles on 
dP. It follows that V^(f has a kernel along V. The kernel of G~^{x) on T^dP is the linear 
span of the normals fir for r G J-'{x). We also denote by H^{p) = VpV9(e'') the Hessian of 
the Kahler potential on the open orbit in p coordinates. By Legendre duality, 

H^{p) = G~\x), p{eP)=x. (57) 

This relation may be extended to V ^ dP. The kernel of the left side is the Lie algebra of 
the isotropy group Gp of any point p G p^^{x). The volume density has the form 



det{G^') = 5^{x)-l[ir{x), {51 



z 



r=l 

for some positive smooth function 6^ [A]. We note that log nr=i ■^^■(^) known in convex 
optimization as the logarithmic barrier function of P. 

2.3. Kahler potential near V. We also need smooth local Kahler potentials in neighbor- 
hoods of points Zq G V. We note that the open orbit Kahler potential fl50l) is well-defined 
near z = 0. Local expressions for the Kahler potential at other points of V essentially amount 
to making an affine transformation of P to transform a given corner of P to 0, and in these 
coordinates the local Kahler potential near any point of V can be expressed in the form fl50|) . 
For instance, on CP^, a Kahler potential valid at 2; = oo is given in the coordinates 10 = ^ 
by log(l -|- |wp). It differs on the open orbit from the canonical Kahler potential log(l -|- \z 
by the term log {zl"^ whose idd is a delta function at 2; = 0, supported on V away from the 
point w = that one is studying. In |So] the reader can find further explicit examples of 
toric Kahler potentials in affine coordinate charts. Hence, in what follows, we will always use 
(!50|) as the local expression of the Kahler potential, without explicitly writing in the affine 
change of variables. 

We will however need to be explicit about the use of slice-orbit coordinates Zj,p" 
in the local expressions of the Kahler potential. The coordinates near zq depend on Tei^o) 
from fl37j) . For each zq & V corresponding to a codimension r face of P, after an affine 
transformation changing the face to x' = 0, we may write the Kahler potential as the 
canonical one in slice-orbit coordinates, F{\z'\'^,e'' ) §2.11 psj) . Since G P, -F is smooth up 
to the boundary face z' = 0. The fact that F is smooth up to the boundary also follows from 
the general fact that a smooth T™-invariant function g G G!^m{M) may be expressed in the 
form g{z) = Fg{p^{z)) where as Fg G C°°(M"*). This is known as the divisibility property 
of T™-invariant smooth functions (cf. [LT] ) . It implies that P is a smooth function of the 
polar coordinates r| near points of P in the sense of fl49p . 

2.4. Almost analytic extensions. In analyzing the Bergman/Szego kernel and the func- 
tions (!2T]) . we make use of the almost analytic extension if{z,w) to M x M of a Kahler 
potential for a Kahler lj; for background on almost analytic extensions, see |BSj MSj] . It is 
defined near the totally real anti-diagonal {z,z) & M x M by 

^,(. + + 5: g^W^^. (59) 
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When if is real analytic on M, the almost analytic extension ip{z,w) is holomorphic in z 
and anti-holomorphic in w and is the unique such function for which '^{z) = ip{z,z). In 
the general C°° case, the almost analytic extension is a smooth function with the right 
side of (l59l) as its C°° Taylor expansion along the anti-diagonal, for which d(f{z,w) = to 
infinite order on the anti-diagonal. It is only defined in a small neighborhood (M x M)s = 
{{z,w) : d{z,w) < 6} of the anti-diagonal in M x M, where d{z,w) refers to the distance 
between z and w with respect to the Kahler metric u. It is well defined up to a smooth 
function vanishing to infinite order on the diagonal; the latter is negligible for our purposes 
(cf. Proposition 1.1 of |BSj| .) 

The analytic continuation (p{z,w) of the Kahler potential was used by Calabi [Caj in the 
analytic case to define a Kahler distance function, known as the 'Calabi diastasis function' 

D{z, w) := (p{z, w) + (p{w, z) - {(^{z) + (p{w)). (60) 

Calabi showed that 

D{z,w) = d{z,wf + 0{d{z,wY), ddlD{z,w%=^ = u. (61) 

One has the same notion in the almost analytic sense. 

The gradient of the almost analytic extension of the Kahler potential in the toric case 
defines the almost analytic extension ^c{z^ w) of the moment map. We are mainly interested 
in the case where w = e^^z lies on the T'^-orbit of z, and by fl53p we have, 

tfic{z,e''z) = VeMz,e''z) = VeFc{e''\z\^), (62) 

where F is defined in (1511) . We sometimes drop the subscript in and fic since there is 
only one interpretation of their extension; but we emphasize that ip{z,e^^z) = Fc(e*^|2p) 
is very different from if{e^^z) = F(|e*^zp) = -F(|2p). For example, the moment map of 
the Bargman-Fock model (C"*, |2:p) is fi{z) = {\zi\'^ , . . . ,\zm\'^), whose analytic extension 
is {ziWi, . . . ,ZmWm)- Similarly that of the Fubini-Study metric on CP*" is (in multi-index 
notation) fiFS,c{z, w) = ^^'^^ . In §2.6l we further illustrate the notation in the basic examples 
of Bargmann-Fock and Fubini-Study models. We also observe that fl62|) continues to hold 
for the Kahler potential F{\z'\'^,e'' ) in slice-orbit coordinates. That is we have, 

t /x(z',e^"/2) = Ve',e"Fc(e^''kT,e^'"+''")|(.',e")=(o,o). (63) 

The complexified moment map is a map 

/ic (M X M)s C™. (64) 

The invariance of /i under the torus action implies that fic{e^^z,e'^^w) = ^c{z,w). The 
following Proposition will clarify the discussion of critical point sets later on (see e.g. Lemma 

Proposition 2.1. For 6 sufficiently small so that fj,c{z,w) is well-defined, we have 

(1) Qfic{z,e''z) = \VeD{z,e''z). 

(2) iic{z, e^^z) = fic{z, z) with {z, e^^ z) e (M x M)s if and only if e''^z = z. 

Proof. The proof of the identity (1) is immediate from the definitions; we only note that 
the diastasis function is a kind of real part, and that the imaginary part originates in the 
factor of i in fl62l) . One can check the factors of i in the Bargmann-Fock model, where 
IJi<c{z,e^^ z) = e*^|2p while D{z,e'^^z) = 2(cos6' — l)\z\'^ + 2i{sm9)\z\'^ (in vector notation) . 
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By (EH), D{z,w) has a strict global minimum aX w = z which is non-degenerate. It is 
therefore isolated for each z. Since its Hessian at w = z is the identify with respect to u, 
the isolating neighborhood has a uniform size as z varies. Thus, there exists a 5 > so that 
fici^y w) = ficiz, z) in (M X M)s if and only ii z = w. This is true both in the real analytic 
case and the almost-analytic case. 

□ 



2.5. Hilbert spaces of holomorphic sections. On the 'quantum level', a toric Kahler 
variety (M, u) induces the sequence of spaces H^{M, L^) of holomorphic sections of powers of 
the holomorphic toric line bundle L with Ci(L) = ^[uo]. The (C*)™ action lifts to H°{M, L^) 
as a holomorphic representation which is unitary on T*". Corresponding to the lattice points 
a G kP, there is a natural basis {s^} (denoted Xa |STZ1] ) of H^{M,L^) given by joint 
eigenfunctions of the (C*)™ action. It is well-known that the joint eigenvalues are precisely 
the lattice points fl kP in the kth dilate of P. On the open orbit Sa{z) = Xa{z)e^ where 
e is a frame and where as above Xa{z) = z°' = z"^ ■ ■ ■ z'^ . Hence, the are referred to as 
monomials. For further background, we refer to [STZlj . A hermitian metric h on L induces 
Hilbert space inner products on H^{M, L'^). 

As is evident from fl?I]) . we will need formulae for the monomials which are valid near P. 
By (go]) and (glD, we have 

Xa.{z) = v,{^)xAz), ze{C*r, (65) 

and by (143|) we then have 

|x.(z)P = |r/^('^f. (66) 

As mentioned above, for simplicity of notation we suppress the transformation F and coordi- 
nates 7], and we will use the 'orbit-slice' coordinates of fH8l) . Thus, we denote the monomials 
cooresponding to lattice points a near a face F by (^z')'^' e^^^^"'^''" where r(a) = {a', a") 
with a" in the coordinate hyperplane corresponding under F to F and with a' in the normal 
space. 



2.6. Examples: Bargmann-Fock and Fubini-Study models . As mentioned above the 
Bargmann-Fock model is the linear model. It plays a fundamental role in this article because 
it provides an approximation for objects on any toric variety on balls of radius and also 
near V. Although it and the Fubini-Study model are elementary examples, we go over them 
because the notation is used frequently later on. 

The Bargmann-Fock models on C"* correspond to choices of a positive definite Hermitian 
matrix H on C'". A toric Bargmann-Fock model is one in which H commutes with the 
standard action, i.e., is a diagonal matrix. We denote its diagonal elements by Hjj. The 
Kahler metric on is thus iddipBF,H{z) where the global Kahler potential is 

m 

V^bf,h{z) = ^Hfj\zj\'^ = F(|2;l|^ . . . , \zm\^), with F{yi, . . .,ym) = ^Hjjyj. 
i=i j 

For simplicity we often only consider the case H = I. Putting {zjl"^ = e''^ and using fl53|l . 
it follows that fiBF,H{zi, . . . , z^) = (i7n|zi|2, . . . , i/^^jz^p) : C"' ^ R^^ as in The 
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symplectic potential Legendre dual to (Pbf,h is given by 



ubf,h{x) = + 2(log/i5^j:^(x),x) = - ^xj + ^Xj\og{j^). (67) 

In this case, Gbf,h is the diagonal matrix with entries ^7-7^-, so detGBF,H = d^^jj'- 
The off-diagonal analytic extension of the Kahler potential in the sense of (l59l) is then 

m. 

^BF,h{z, w) = HfjZjWj = F{ziWi, . . . , ZmWm) 

i=i 

and in particular, 

m 

^BF,H{^,e''z) = Y,Hf,e'''\^3\' = He'''\zi\', • • .,e''\zm\'). 

i=i 

Henceforth we often write the the right side in the multi-index notation Fc{e^^\z\'^). We 
observe, as claimed in fl^ . that Ve-FBF,c(e*^kP)|e=o = ifJ'BFi^). 

Quantization of the Bargmann-Fock model with H = I produces the Bargmann-Fock 
(Hilbert) space 

H'(C'", (27r)-'"A;"^e-^l"l'dz A dz) 

of entire functions which are relative to the weight e~^^^^^^^. It is infinite dimensional and a 
basis is given by the monomials z"" where a G W^dl/^. In §3.0.1l we compute their norms. 
For H ^ I one uses the volume form e'''^^^'^^ {idd{Hz, z))"^ /m\ = e~'^^^^'^^(det H)dz A dz. 

Toric Fubini-Study metrics provide compact models which are similar to Bargmann-Fock 
models. In a local analysis we always use the latter. A Fubini-Study metric on CP"* is 
determined by a positive Hermitian form H on C"*"*"^ and a toric Fubini-Study metric is a 
diagonal one ^JIq In the affine chart Zq ^ (e.g.) a local Fubini-Study Kahler 

potential is lpfs,h{zi, . . . , Zm) = log(l + hj-^\zj\^) where hfj = jf^. This is a valid Kahler 
potential near z = but of course has logarithmic singularities on the hyperplane at infinity. 
The almost analytic extension of the Fubini-Study Kahler potential is given in the affine 
chart by log(l -|- hfjZjWj). Thus ( j62l) asserts that 
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Velog(l + ^/i 



l2\ 



Quantization produces the Hilbert spaces H°{CF"^,0{k)), where 0{k) CP"* is the 
kth power of the hyperplane section bundle. Sections lift to homogeneous holomorphic 
polynomials on C""*"^, and correspond to lattice points in /cS where E is the unit simplex in 

R"\ 

2.7. Linearization of the Monge- Ampere equation. It is known that the Legendre 
transform linearizes the Monge-Ampere geodesic equation. Since it is important for this 
article, we present a simple proof that does not seem to exist in the literature. 

Proposition 2.2. Let Mp be a toric variety. Then under the Legendre transform ip — > u^, 
the complex Monge-Ampere equation on Tix™ linearizes to the equation u" = 0. Hence the 
Legendre transform of a geodesic ipt has the form Ut = uq + t{ui — uq). 
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Proof. It suffices to show that the energy functional 

E = [ [ ^Idii^^dt (68) 

Jo JM 

is Euchdean on paths of symplectic potentials. For each t let us pushforward the integral 
Jj^jip'fdfi^ under the moment map /i^^. The integrand is by assumption invariant under the 
real torus action, so the pushforward is a diffeomorphism on the real points. The volume 
measure d/i^j pushes forward to dx. The function dtipt{p) pushes forward to the function 
ipt{x) = ^t{Px,t) where PipXPx,t) = x. By fl55l) . the symplectic potential at time t is 

ut{x) = {x,px,t) - ^t{Px,t)- 

We note that 

iit = {x, dtPx,t) - 'Pt{Px,t) - {'^pVt{Px,t), dtpx,t)- (69) 
The outer terms cancel, and thus, our integral is just 

/ / \ut\'^dx. 
Jo Jp 

Clearly the Euler-Lagrange equations are linear. 

□ 



3. The functions Vhk and Q^k 



We now introduce the key players in the analysis, the norming constants Qhk{a) ( !20l) and 
the dual constants Vh^{a) of ( l22l) . The duality is given in the following: 



Proposition 3.1. We have: 



Proof. By it follows that 



Kip-.'i^ml = lx«(/i.-^(^))|V^-('^r(f)) = e^'^-(f). (70) 



Corollary 3.2. 



7^fc(^,a) 



Proof. We need to show that 



□ 



(71) 



By Proposition 13. ![ the left side of fl7T|) equals 
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By (CD]), the left side of ([71]) equals 



But Ut{x) + (1 — t)no(x) + tui{x) = on a toric variety, and this gives the stated equality. 

□ 

Further, we relate the full Vhk{a,z) to the Szego kernel. The Szego (or Bergman) kernels 
of a positive Hermitian line bundle (L, h) — >■ (M, u) over a Kahler manifold are the kernels 
of the orthogonal projections U^k : L'^{M, L^) H^{M, L^) onto the spaces of holomorphic 
sections with respect to the inner product Hilhkih) ([6]). Thus, we have 



Jm ml 

where the ■ denotes the /i- hermitian inner product at w. Let ej;, be a local holomorphic frame 
for L — ^> M over an open set f/ C M of full measure, and let {sj = fjcf^ : j = 1, . . . , dk} he 
an orthonormal basis for H^{M,L^) with = dim if°(M, L'^). Then the Szego kernel can 
be written in the form 



U,k{z,w) ■.= F,k{z,w)ef{z)(^ef{w), (73) 

where 

dk 

F,k{z,w) = Y,fji^)W)- (74) 
i=i 

Since the Szego kernel is a section of the bundle (L'^) ® {L'^)* — > M x M, it often simplifies 
the analysis to lift it to a scalar kernel tl^k {x, y) on the associated unit circle bundle X M 
of {L,h). Here, X = dDJ^ is the boundary of the unit disc bundle with respect to in 
the dual line bundle L*. We use local product coordinates x = {z,t) E M x on X where 
X = e^^\\eL{z)\\he*];^{z) G X. To avoid confusing the action on X with the action on M 
we use e** for the former and e*^ (multi-index notation) for the latter. We note that the T"* 
action lifts to X and combines with the 5*^ action to produce a {S^)^~^^ action. We refer to 
[Zell ISZ| IZe2] for background and for more on lifting the Szego kernel of a toric variety. 

The equivariant lift of a section s = fef' G H'^{M, L^) is given explicitly by 

s{z^) = e^^%efUf{z) = e4-i^WH/(^) . (75) 
The Szego kernel thus lifts to X x X as the scalar kernel 

tik{z,t-w,t') = e*^[-^'^(^)-^^(")+'(*-*')]Ffc(2,w) . (76) 
Since it is S^- equivariant we often put t = t' = 0. 
Proposition 3.3. We have 

Vhk{a,z) = (27r)-™ [ 11^^(6''^ z,0; z, 0)6-'^"''^^ d9. 
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Proof. We recall that Xa{z) = 2° is the local representative of Sa in the open orbit with 
respect to an i 

basis, we have 



respect to an invariant frame. Since {—7^==} is the local expression of an orthonormal 



Qhk(a) 



hence 



Qhk(a) 
It follows that 

\Xa{z)\'^e-^'P^'^e'^''^^^ 



aGfcPn: 

Integrating against e~*^"'^^ sifts out the a term. 



□ 



Corollary 3.4. We have 



= (27r)-- / n,.(e^V.-'(^),0;/i,-^(^),0)e-<"'^>d^. (77) 



3.0.1. Bargmann-Fock model. As discussed in §2.6[ the Hilbert space in this model has the 
orthogonal basis z" with a G M™nZ™. The Bargmann-Fock norming constants when H = I 
are given by 

g,|^(a) = fc-H— «!, {a\:=a,\---aj.) 

It follows that an orthonormal basis of holomorphic monomials is given by {/c' 'J -^}- 
We therefore have 

— ^ = A;l°l+-^e-'l^l , (78) 



and in particular, 

(a) = A;™e-H^, (79) 

where a" = 1 when a = 0. Here, we use that ubf{^) = f logf ~ so that e'^"^^'-^^ = 
g-|a|fc_^ and that Q^fe^(a) = k~'^~^"^a\. We observe that Vfik^^{a) depends on k only 
through the factor A;™". 

Precisely the same formula holds if we replace / by a positive diagonal H with elements 
Hjj. By a change of variables, 2/i|^^(a) = ^]LiH ^-""^ Qf^k^^{a) , and also by (1^ ubf,h{x) = 
ubf{x) + J2j ^j^^sHjj. Hence, by Proposition 13. ![ 

V,,k {a)=Vf,k {a)U"l.H~-"'e^^''^^°^^f^ =Vf,k (a). 
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3.0.2. CP™'. In the Fubini-Study model, a basis of H^{CF"^,0{k)) is given by monomials 
with a G fcS (see §2.61) . and the norming constants are given by 

Q 



ihk [a) 



«1, 



0) 



Recall that multinomial coefficients are defined for ai 



a.m < by 



Q!l, 



ar. 



where as above, |a| = ai + 



We further have |sa(2;)|^fc 

"■FS 



.ai2g-fciog{i+|^|2) therefore. 



and since 



we have 



k\ 



\a\ 

T 



,fc-|a| 



'1- 



\a\ 

T 



,fc-|a| 



ai! ■ ■ ■ ctmK^ ~ l"^!)' ^ 

4. SZEGO KERNEL OF A TORIC VARIETY 

We will use Proposition 13. 3l to reduce the joint asymptotics of Vh^a, z) in (/c, a) to asymp- 
totics of the Bergman-Szego kernel off the diagonal. We now review some general facts 
about diagonal and off-diagonal expansions of these kernels, for which complete details can 
be found in |SZ] , and we also consider some special properties of toric Bergman-Szego kernels 
which are very convenient for calculations; to some extent they derive from [STZlj . but the 
latter only considered Szego kernels for powers of Bergman metrics. 

The Szego kernels ti^ki^x, y) are the Fourier coefficients of the total Szego projector n/i(x, y) : 
— > li?{X), where li?{X) is the Hardy space of boundary values of holomorphic func- 
tions on D* (the kernel of db in L'^{X)). Thus, 

1 [^'^ 

Yihu (x, y) = — e^'^'lih{e''x, y) dt. 
27r Jo 

The properties we need of Yif^k{x,y) are based on the Boutet de Monvel-Sjostrand construc- 
tion of an oscillatory integral parametrix for the Szego kernel ( |BSj] ): 



Ii{x,y) = S{x,y) + E{x,y) 



^1) 



with S{x, y) = e*^^(^'^)s(x, y, X)dX , E{x, y) G C°°(X x X) . 
The phase function ip is of positive type and is given in the local coordinates above by 

i){z, t; w, t') = -\l- e^(^''")-|(^W+^{'^))e*(*-*') . 
i L 

Here, ip{z,w) is the almost analytic extension of the local Kahler potential with respect 
to the frame, i.e., h = e"'^'-^-'; see fl59l) for the notion of almost analytic extension. The 



2) 
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amplitude s(^z,t;w,t' , Xj is a semi-classical amplitude as in |BSj] (Theorem 1.5), i.e., it 
admits a polyhomogeneous expansion s ~ J^JLq \"^~''sj{x, y) G S^{X x X x IR+). 

The phase '^(2;, t; w , t') is the generating function for the graph of the identity map along 
the symplectic cone E C T*X defined by S = {(x, raj.) : r > 0} where is the Chern 
connection one form. Hence the singularity of II (x, y) only occurs on the diagonal and 
the symbol s is understood to be supported in a small neighborhood (M x M)^ of the anti- 
diagonal. It will be useful to make the cutoff explicit by introducing a smooth cutoff function 
x{d{z^ w)) where x is a smooth even function on M and d{z, w) denotes the distance between 
z, w in the base Kahler metric. 

As above, we denote the k-th Fourier coefficient of these operators relative to the 
action by Hj^^k = Si-^k + E^k. Since E is smooth, we have Ef^k{x, y) = 0{k~°°), where 0{k~°°) 
denotes a quantity which is uniformly 0{k~"') on X x X for all positive n. Hence Ef^k{z, w) is 
negligible for all the calculations and estimates of this article, and further it is only necessary 
to use a finite number of terms of the symbol s. For simplicity of notation, we will use the 
entire symbol. 

It follows that (with x = (z,t),y = (w, 0) and with x{d{z,w)) as above ), 
flh^ix,y) = Shkix,y) + Oik-°°) 

(83) 

= ir IS" e^^'(-*+^^(^'*'"''°)) x{d{z, w)) s{z, t; w, 0, k\)dtd\ + 0{k~°°) 

The integral is a damped complex oscillatory integral since 061 p implies that 

Qi:{x,y)>Cd{x,yY, {x,y e X), (84) 

for [x, y) sufficiently close to the diagonal, as one sees by Taylor expanding the phase around 
the diagonal (cf. |BSj| , Corollary 1.3). It follows from fl55]l and from flMl) that the Szego 
kernel Ilhk{z,w) on M is 'Gaussian' in small balls d{z,w) < i.e., 

\Ilhk{z,^;w,ip')\ <Ck"'e-^'^^'^'"^' + 0{k-°°), {when d{z,w) <^-^), (85) 

and on the complement d{z, w) > it is rapidly decaying. This rapid decay can be 
improved to long range (sub-Gaussian) exponential decay off the diagonal given by the 
global Agmon estimates, 

\nh^{z,ip;w,ip')\ < CA;™e"^'^("'"'^ (86) 

We refer to [Chrl [L] for background and references. 

It is helpful to eliminate the integrals in fl83l) by complex stationary phase. Expressed in 
a local frame and local coordinates on M, the result is 

Proposition 4.1. Let {L,h) be a C°° positive hermitian line bundle, and let h = e""^ 
in a local frame. Then in this frame, there exists a semi-classical amplitude Ak{z,w) ~ 
k'^ao{z, w) + k"^~^ai{z, w) + ■ ■ ■ in the parameter k~^ such that, 

n,.(z, 0; w, 0) = e^('^(^''")-5('^W+'^(-)))xfc(d(z, w)) Ak{z, w) + 0{k-^), 

where as above, Xk{d{z,w)) = x{^-^d{z,w)) is a cutoff to - neighborhood of the diagonal. 
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Proof. This follows from the scaling asymptotics of [SZ] or from Theorem 3.5 of [?]. We refer 
there for a detailed proof of the scaling asymptotics and only sketch a somewhat intuitive 
proof. 

The integral (!83|) is a complex oscillatory integral with a positive complex phase. With 
no loss of generality we may set (f' = 0. Taking the A-derivative gives one critical point 
equation 

X _ g'/'(2."')-|('/'(2)+'^(«'))g«i _ g 

and the critical point equation in t implies that A = 1. The A-critical point equation can only 
be satisfied for complex t with imaginary part equal to the negative of the 'Calabi diastasis 
function' (1^ . i.e., 

St = D{z, w), 

and with real part equal to —Q(f{z,w). To obtain asymptotics, we therefore have to deform 
the integral over to the circle \(\ = e"^^^'"'). Since d{z,w) < C-^^ by assumption, the 

deformed contour is a slightly re-scaled circle by the amount in the complete proofs, 
the contour is held fixed and the integrand is rescaled as in |SZ] . The contour deformation 
is possible modulo an error 0{k~^'^) of arbitrarily rapid polynomial decay because the inte- 
grand may be replaced by the parametrix (up to any order in A) which has a holomorphic 
dependence on the C* action on L*, hence in e*^ to a neighborhood of in C This is 
immediately visible in the phase and with more work is visible in the amplitude (this is the 
only incompleteness in the proof; the statement can be derived from |SZj and also [Chrj ). 
We need to use a cutoff to a neighborhood of the diagonal of M x M, but it may be chosen 
to be independent of 9. 

By deforming the circle of integration from the unit circle to |C| = e^^^''^-' and then changing 
variables t t + iD{z, w) to bring it back to the unit circle, we obtain 

poo p27T 

Ilhk{x,y)^k / / e''=(-*-'^("''")-^'^("'*+'^("''")''"'°))s(^,t+2D(z,w^ mod fc"^ 

Jo Jo 

The new critical point equations state that A = 1 and that q^'=^'p{z,w) ^tt _ calcu- 
lation shows that ip = Q o\i the critical set so the phase factor on the critical set equals 
^ip{z,w)-\{ip{z)+ip{w)) ^ rjj-^g Hessian of the phase on the critical set is ( ? ^ | as in the di- 

agonal case and the rest of the calculation proceeds as in [Zelj . (As mentioned above, a 
complete proof is contained in [SZ] ). 

□ 

4.1. Toric Bergman-Szego kernels. In the toric case, we may simplify the expression for 
the Szego kernels in Proposition 14. II using the almost analytic extension (cf. §2.41) : ( l59l) ) of 
the Kahler potential ip{z.,w) to M x M, which has the form 

F<c{z ■ w) = the almost analytic extension of -F(|2;|^) to M x M. (88) 

The almost analytic extension will be illustrated in some analytic examples below, where it 
is the analytic continuation. 
Thus, we have: 
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Proposition 4.2. For any hermitian toric positive line bundle over a toric variety, the 
Szego kernel for the metrics have the asymptotic expansions in a local frame on M , 

where Ak{z^ w) ~ fc™ [cioiz^ w) + + ■ ■ ■ j is a semi- classical symbol of order m. 

As an example, the Bargmann-Fock(-Heisenberg) Szego kernel with k = 1 and H = I is 
given (up to a constant Cm depending only on the dimension) by 

cel 

The higher Szego kernels are Heisenberg dilates of this kernel: 

fl^,^(a;,y) = -l-fcV'=(*-^)e'=(«-*'-^l^l'-^l''l'), (89) 

where x = {C,t) , y = [t], s) . In this case, the almost analytic extension is analytic and 
Fbf,c{z,w) = z -w. 

A second example is the Fubini-Study Szego kernel on 0{k), which lifts to 5^™"^ x S*^™"^ 

as 

U.k (x,y)=y ^—-x-^y=- Tv-{x,y) . (90) 

Recalling that x = e^^TT^^n in a local frame e over an affine chart, the Szego kernel has the 

||e(2)|| ' fe' 

local form on C™ x C™ of 

n,.^(^,0;^,0) = ^^^^e ^v^vOT. (91) 
Thus, Ffs,c{z, w) = log(l + z ■ w). 

4.2. Asymptotics of derivatives of toric Bergman/Szego kernels . One of the key in- 
gredients in of Theorem 1 1.1 1 is the asymptotics of derivatives of the contracted Bergman/Szego 
kernel 

n^. (z, z) = F,. (z, z) I (z) 1 1^. = (z, 0; z, 0) (92) 

in {t,z). (The notation is slightly ambiguous since in fl73l) it is used for the un-contracted 
kernel, but it is standard and we hope no confusion will arise since one is scalar-valued and 
the other is not.) These derivatives allow us to make simple comparisions to derivatives 
of ipk{t,z). Since we ultimately interested in norms we need asymptotics of derivatives 
with respect to non-vanishing vector fields. We can use the vector fields ^ away from V 

and the vector fields 7^ near T). The calculations are very similar, but we carry them both 
out in some detail here. Later we will tend to suppress the calculations with ^ to avoid 
duplication; the reader can check in this section that the calculations and estimates are valid. 

Only the leading coefficient and the order of asymptotics are relevant. The undifferentiated 
diagonal asymptotics are of the following form: for any h G P{M,u), 

dk 

U^,{z, z) = J2\\si{z)\\l = aok^ + a,{z)k^-^ + a2{z)k^-' + ... (93) 

i=0 
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where oq is constant and as above dk + 1 = dim^r°(M, L''). 

We first consider derivatives with respect to p. Calculating p derivatives of n^fe(e^/^, e^/^) 
is equivalent to calculating ^-derivatives of Il,^k{e^^z, z). Using flB^ we have 

n,.(e^^.)= — 



The results are globally valid but are not useful near T) since on each stratum some of the 
vector fields generating the (C*)"^ action vanish. 

In the following, we use the tensor product notation (f — /Ut(e^''^))^-^ for (^— /ii(e''/^)j)(^ — 

Proposition 4.3. For i,j = l,...,mwe have, 

(2) ri;;^ (-E...Pn.4|logQ..(a)) ^^^^^^Jf^) - ^1^* = 0(A:-); 

(4) n;;^ (feE...Pnz4f - /i.(e^/^)).(|log Q..^(c.)) -^°t,;(irO " ^gt = ^(fe"^)- 
Proof. To prove (1), we differentiate and use (l53 l) -( l62l) and ( l93i) to obtain 

To prove (2) we differentiate 

log^,.(e''/^e''/2)=log J] ^ , , 

with respect t to produce the left side. Since the leading coefficient of (l93l) is independent 
of t, the t derivative has the order of magnitude of the right side of (2). 
To prove (3), we take a second derivative of (1) in p (or 6) to get 



Then (3) follows from (l93ll and the fact that V/it(e^''^)) = V^(^. Similar calculations show 
(4). 

□ 

In our applications, we actually need asymptotics of logarithmic derivatives. They follow 
in a straightforward way from Proposition 14.31 using that Ilhk{z,z) ~ k"^. We record the 
results for future reference. 

Proposition 4.4. We have: 
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iV loeV ^- $—-0(^] 



|5c(2)|2, JZadtlogi ) 



VI) 



Proposition 4.5. We have: 

\S f^lP / \ ~^ 



(q) / I dpidp 



(3) 



e,fc(«) 



|Sc«(2)| 



1 ^2 , '"^'~"hk 

k 5*2 Z^aGfcPnZ™ Q . (a) 

"t 

E.., (d? log (a log i°g(^^)) e^^'g) , , 

Eq S fe'(Q) 

Finally, we consider the analogous derivatives with respect to the radial coordinates rj 
near V. We assume z is close to the component of V given in local slice orbit coordinates 
by 2;' = and let r' = (vj)^^^ denote polar coordinates in this slice as discussed in ^ The 
Szego kernel then has the form 

r'^°'J ^{p" ^-kFt{rl...,rj,e''P+\...,ePm) 

Uj,,(z,z)= y -^^^ — — . (94) 

The coefficients of the expansion (p3l) are smooth functions of r| and the expansion may be 
differentiated any number of times. 

The behavior of n^fc(^,^) for z & T) has the new aspect that many of the terms vanish. 
The extreme case is where z is a fixed point. We choose the slice coordinates so that it has 
coordinates 2 = 0. We observe that only the term with a = in (IMll is non-zero, and the 
ath term vanishes to order 

Since S- = —jr- where both are defined, the calculations above are only modified by the 

drj rj dpj ' J 

presence of new factors of ^ in each space derivative. Since we are applying the derivative to 

functions of r|, it is clear that the apparent poles will be cancelled. Indeed, the Vj derivative 
removes any lattice point a with vanishing aj component. Comparing these derivatives with 
derivatives of flMl) gives the following: 



Proposition 4.6. For n = 1, . . . ,p, we have: 
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T,P to" a") -kFt[r^,...,rl,e^V+^ ....,ePm) 



In effect, the exponent a is taken to a — (0, . . . , 1„, . . . ) in the sum or removed if = 0, 
where (0, . . . , 1„, . . . ) is the lattice point with only a 1 in the nth coordinate. There are 
similar formulae for the second derivatives 73-^73—, tt-ju^ a '^a ■ The only important point to 
check is that the modification changing a to a — (0, ) does not affect the proofs in 



5. Localization of Sums: Proof of Lemma [TT2] 
The following Proposition immediately implies Lemma \1.2\ 
Proposition 5.1. Given {t,z), and for any 6,C > 0, there exists C > such that 

- VAa,z) = Oik-''), if 1^ - ^^t{z)\ > C'k~h+^. 



Proof. The proof is based on integration by parts. All of the essential issues occur in the 
Bargmann-Fock model, so we first illustrate with that case. 

5.L Bargmann-Fock case. To analyze the decay of Vhk^^{a, z) as a function of lattice 
points a, it seems simplest to use the integral formula (suppressing the factor fc™' and nor- 
malizing the volume of T™ to equal one), 

(95) 

We observe that the rightmost expression in fl95l) is e times a complex oscillatory 
integral with phase 

We observe that (consistent with Proposition 12. ip . 

V,$,.(^) = z(k|V^ - f ) = ^ = 1^1' = I- 



' k 

Further, we claim that 



k' 



Indeed, the function 



on T™ has a strict global minimum at 6' = as long as 7^ 0, 7^ for all j. It still has 
a global minimum without these restrictions, but the minimum is no longer strict. We note 
that this discussion of global minima is possible only because the Kahler potential admits a 
global analytic continuation in {z,w)] in general, one can only analyze critical points near 
the diagonal. 



BERGMAN METRICS AND GEODESICS 



33 



We integrate by parts with the operator 



1 1 



k |Ve$^,||2 



i.e., we apply its transpose 



to the amphtude. The second (divergence) term is —1 times 



(99) 



k |V$^-P k |V$^-|4 

I ^ k ' ' ' fc ' 

We will need to take into account the /c-dependence of the coefficients, and therefore 
introduce some standard spaces of semi-classical symbols. We denote by 5'^(T™') the class 
of smooth functions ak{0) on T™ x N satisfying 

sup I^X-WI < (100) 

Here we use multi-index notation Dg = YYjLiijio')'^^ ■ Thus, each D^. derivative gives rise to 
an extra order of k^ in estimates of a^. We note that products of symbols satisfy 

S^' X c ^5"^+"^ (101) 
We now claim that (with 6 the same as in the statement of the Proposition), 

(2) (I99D lies in S^^^. 

2 

In (2), we note the pre-factor i. To prove the claim, we ffist observe that the sup norm 

estimates are correct by fl96l) and from the fact that |^^"'^ | is a unit vector. We further 

consider derivatives of (l)-(2). Each 6 derivative essentially introduces one more factor 
of k\'Vg^z,f\ and hence raises the order by /c2~^. This continues to be true for iterated 
derivatives, proving the claim. 
Now we observe that 

C : ^ Slzf. (102) 

2 2 

Indeed, the ffist term of is the composition of (i) Vg, which raises the order by | — 5, 

(ii) multiplication by an element of 'S'i_^ which again raises the order by | — 5 (iii) times | 

which lowers the order by 1. The second term is a multiplication by ^ times an element of 
Sl^'^f and thus also lowers the order by 26. 

It follows that each partial integration by C introduces decay of k~'^^, hence for any M > 0, 



m = e-^l^l'(27r)— /^„eKl^l'^")-<f'^))((/:*)^^l)d^ 



= U(fc— 

in this region. 
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5.2. General case. We now generalize this argument from the model case to the general 
one. With no loss of generality we may choose coordinates so that z lies in a fixed compact 
subset of C^, where the open orbit is identified with (C*)™'. In the open orbit we continue to 
write \z\^ = e''. The first step is to obtain a useful oscillatory integral formula for Vh^iot, z). 
By Proposition 13.31 and Proposition 14.21 we have 

Vhk{a,z) = (27r)-™ f e'=(^c(e»'^I^P)-nkP)) ^{^d{z , e''^ z)) Ak{z, e'^ z,0)e'^'''^^ dO +0(A;^°°). 

(103) 

The phase is given by 

<l>.,|(^) = Fc(e*^|^n - F{\z\^) - 6) (104) 

where as above, Fc(e*^|2;p) is the almost analytic continuation of the Kahler potential F(|2p) 
to M X M. By (IHID and ([6l]), it satisfies 

^{Fc{e'^\z\'') - Fd^H)) < -Cd{z, e'^zf, (for some C > 0). (105) 

Hence, the integrand fll03p is rapidly decaying on the set of 6 where d{z,e^^z)'^ > (7^^ 
(see also (186|1 ). and we may replace xi^iz, e^^ z)) by xik^~^'d{z,e^^z)) G since the 

2 " 

contribution from 1 — x(^^ d{z,e'^^z)) is rapidly decaying. Here, 6' is an arbitrarily small 
constant and we may choose it so that 6' < 6 in the Proposition. (We did not use such cutoffs 
in the Bargmann-Fock case since the real analytic potential had a global analytic extension 
with obvious properties, but as in §2.41 it is necessary for almost analytic extensions). 

The set d{z, e^^z) < depends strongly on the position of z relative to P, or equivalently 

on the position of ^h{z) relative to dP. For instance, if 2; is a fixed point then d{z^ e^^z) = 
for all 6. However, we will not need to analyze these sets until the next section. 

We now generalize the integration by parts argument. Our goal is to prove that Vf^k (a, z) = 

0{k~'-') if If — l^t{z)\ > Ck~^~^^. Now, the gradient in 9 of the phase of fll03|] is given by 

Ve^z, = VeFc{e^'\z\') - = z{fic{z, e^'z) - ^), (106) 

where ^c{z,e^^z) is the almost analytic extension of the moment map (see §2.4p . The fol- 
lowing Lemma is obvious, but we display it to highlight the relations between the small 
parameters 6 of the Proposition and 6' in our choice of cutoffs. 

Lemma 5.2. // |f - fit{z)\ > Ck-i+^, and if d{z,e'^z) < Ck-^+^' with S' < S , then 
\{fi{z,e''^z) - f)| > C'k--2+\ 
Proof. By Proposition 12.11 

l^^^,^e''z)-f)\' = \{sfi^{z,e''z)-^)\' + \lVoD{z,e''z)\' 

> \(^(z)-^)\^ + 0{d{e^'z,z)). 

□ 

It follows that under the assumption |^ — fit{z)\ > Ck~^~^^ of the Proposition, we may 
integrate by parts with the operator 

C = i|Ve$.,|r' Ve$,,| ■ V, . (107) 
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The transpose has the same form ( l98l) as for the Bargmann-Fock example, the only 
significant change being that it is now applied to a non-constant amplitude Ak and to the 
cutoff x{k^~^ d{z,e'^^z)) E S'i_,, as well as to its own coefficients. Differentiations of Ak 

2 ° 

do preserve the orders of terms; the only significant change in the symbol analysis in the 
Bargmann-Fock case is that differentiations of x{k^~^ d{z, e^^z)) bring only improvements of 
order rather than k~^. However, the order still decreases by at least 26' on each partial 
integration, and therefore repeated integration by parts again gives the estimate 

\Vh.{a,z)\ = O (^{k''')'' j e'=(«^(^"l^l')-^(l^l')de^ = O {{k~''Y') . (108) 

□ 

Remark: It is natural to use integration by parts in this estimate since the decay in ^t{z) — f 
must use the imaginary part of the phase and is not a matter of being far from the center of 
the Gaussian. 

5.3. Further details on the phase. For future reference (see Lemma 16.21) . we Taylor 
expand the phase (11041) in the 9 variable to obtain 

<I>,,.(^) = ^{^^{z) -^,0) + (if.^,^) + i?3(A:,e V'(^)), (109) 

where R3 = 0{\ef). Here, = V^F{fi-^{^)) denotes the Hessian of at f (see ([57]) of 
§2.2p Indeed, we have 

Fc{e^'\z\') - F{\z\') = iFc{e^''\z\')dt 

= J^{VeF{e''<'\z\'),te)dt 

= {V,F{e^)), {t9)) + J^{t - l)V%F{e^''+P)){te)y2dt (HO) 

= 2(M^),^) + V2(F(eO)(^^)2 + i?3(fc,e>-i^"^ 



= z(/i(^), 6) + {H^e, e) + R;{k, e, /i-^f )), 

in the notation ( 157|) . where = V^pF{\z\'^) and where 



Rs{k,e,p):= / it-iy{V;iF{e'''^^)),{teY/3\)dt. (Ill) 

^0 

6. Proof of Regularity Lemma [TT^ and joint asymptotics of Vhk{a) 
The first statement that IZooit.x) is up to the boundary follows from (l58l) . 



"^OO (^) •^) 



, 1/2 



(112) 



(5^(j(x)i-'5^-^(x)* 
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where the functions 5^, are positive, bounded below by strictly positive constants and C°° 
up to dP. 

We now consider the asymptotics of TZk{t, a). We determine the asymptotics of the ratio 
by first determining the asymptotics of the factors of the ratio. We could use either the 
expression (!30!) in terms of norming constants Q\{(y) for the dual expression in terms of 
Vhk{a) in Corollary 13. 2[ Each approach has its advantages and each seems of interest in the 
geometry of toric varieties, but for the sake of simplicity we only consider Vh^ (a) here. In 
[SoZlj we take the opposite approach of focusing on the norming constants. The advantage 
of using Vhk{a) is that it may be represented by a smooth complex oscillatory integral up 
to the boundary, while Q^(a) are singular oscillatory integrals over P. A disadvantage of 
Vh.k{a) is that it does not extend to a smooth function on P and has singularities on dP. 

The asymptotics of (a) are straightforward applications of steepest descent in compact 
subsets of M\D but become non-uniform at T). To gain insight into the general problem we 
again consider first the Bargmann-Fock model, where by fl79l) we have 

V^k {a) = /t'^e-l"!^ = (27r)^™A;™ / e'"'^''"-^-"^''^^ dO. (113) 

As observed before, the factors of k cancel so 'asymptotics' means asymptotics as a ^ oo. 
This indicates that we do not have asymptotics when a ranges over a bounded set, or 
equivalently when ^ is j-close to a corner. On the other hand, steepest descent asymptotics 
applies in a coordinate aj as long as aj oo. Our aim in general is to obtain steepest 
descent asymptotics of Vhk{a) in directions far from facets and Bargmann-Fock asymptotics 
in directions near a facet. 

6.1. Asymptotics of Vhk{a). The analysis of V^k^a) is closely related to the analysis of 
Vhk{<y, z) in §5.21 and in a sense is a continuation of it. But the arguments are now more than 
integrations-by-parts. We obtain the asymptotics of V^k {a) from the integral representation 
analogous to fll03p (see also Proposition 14.21 and Corollary 13.41) . Modulo rapidly decaying 
functions in /c, we have (in the notation of Proposition 14. 2p . 



V^k{a) ~ {27iy-S^„^e-^^^-^^''^nH-k))-n,lHV^^ 

(114) 

This largely reduces the asymptotic calculation of Vhk{a) to facts about the off-diagonal 
asymptotics of the Szego kernel (cf. Proposition 14. 2p . 

The integral f lll4p is the oscillatory integral (11031) but with z = /i"^(^). Hence, as in 
(11040 . its phase is 

^.{6) = Fc(eV^(^)) - Fif^-'i^)) - ^(p (115) 
As in ( IH^ and (llOSp (but with i included in as part of the phase) , 

m^{e) < -Cd{fi'\j),e''fi~\j))\ (for some C > 0). (116) 
Specializing ( I106P to our z = yU~^(|), we get 

V.$.(^) = V,Fc(eV^(f )) - = ^(/^c(/i-^(f ),eV^(f )) - f )• (117) 
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By Proposition 12. 11 the complex phase has a critical point at values of 9 such that d{z, e^^z) < 
5, and e*^/i~^(^) = For ^ ^ dP, the only critical point is therefore 6' = 0. The phase 

then equals zero, and hence at the critical point the real part of the phase is at its maximum 
of zero. 

For J ^ dP, the critical point 6' = is non-degenerate. Specializing fllOQp to z = 
we have 

^c(e>^^(f))-F(/.-^(f)) = /o'|i^c(e^*V.-^(f))cit 

(118) 

where Rs is defined in f ill II) . Hence, 

(^) = (H^e, 9) + R;{9, k, (119) 

and finally, 

V,.{a) ~ {2nr"^J^^e-'^''^'''^e>'^-^('''''^-'(^^^ . (120) 

Non-degeneracy of the phase is the statement that is a non-degenerate symmetric matrix, 
and this follows from strict convexity of the Kahler potential or symplectic potential, see 
( 157|) . But as discussed in §2.21 the H<^ has a kernel when | G dP. Hence the stationary 
phase expansion is non-uniform for f G P and is not possible when ^ G dP. This accounts 
for the fact that we need to break up the analysis into several cases, and that we cannot rely 
on the complex stationary phase method for all of them. 
Specializing fl85|) and fl86|) . we have 

I u^,^e^'fi^\j),fi-\j))\ < Cfc-e-^'='^(^)'^"t))V0(e-^^'^(^''="^)). (121) 

rh K 

Hence, the integrand of (fTT4D is negligible off the set of ^ where c/(/i-i(f ), e^V^Hf )) < C*^- 
We now observe that for d{z, e^^ z) < C^, 

d{e''^z,zf ~ ^(1 -cos%)£j(/i(^)), (122) 
i 

where we sum over j such that |£j(/i(z))| << 1 (we will make this precise in Definition 16. ip . 
In particular, 

d{e''^^l\l).^^l\l)? - - cos^,)£,(^). (123) 

3 

Indeed, both in small balls in the interior and near the boundary, the calculation is universal 
and hence is accurately reflected in the Bargmann-Fock model with all Hj = 1, where the 
distance squared equals 

mm m 

J2 V^'^3 - ^jf = - cos^j) = 2 J]£j(/i(z))(l - cos^j). (124) 

j=i j=i j=i 

This motivates the following terminology: 

Definition: Let < 5^ << 1 . We say: 

• X G P is (5fc-close to (resp. b^-iw. from) the facet Fj = {ij = 0} if ij{x) < 6k (resp. 
>Sk). 
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• X is a (5fc-interior point if it is 5/c-far from all facets. 

There are m possible cases according to the number of facets to which x is 6k-close. Of 
course, x can be (5fc-close to at most m facets, in which case it is (5fc-close to the corner defined 
by the intersection of these facets. We thus define 

J^s,{x) = {r:\ir{x)\<6k}. (125) 

We also let 

5f(x) = #^4(x) (126) 
denote the number of ^^-close facets to x. Dual to the sets J^s^ above are the sets 

^F,,,...,F,, = {x : J^5,(x) = {zi,...,v}}. (127) 

The asymptotics of Vhk{a) depend to the leading order on the determinant of the inverse 
of the Hessian of the phase of (11141) (see also (11031) ) at ^ = 0. This Hessian is the same as the 
Hessian of the Kahler potential discussed in §2.2[ and we recall that its inverse is the Hessian 
G of the symplectic potential. Hence, the asymptotics are in terms of the determinant of G, 
which has first order poles on dP. This indicates that the asymptotics are not uniform up 
to dP. We saw this as well in the explicit example of the Bargmann-Fock case. We define 

g^a{^)= [w- n ^^■(^) ' (128) 

where the functions 5^ are defined in §2.2[ When x is (5fc-far from all facets, then Q^{x) = 
det (cf. [58]). We also define V^k^^^Mj^x)) to be the unique real analytic extension of (j79l) 
to all X G [0, oo). We then consider Bargmann-Fock type functions of type fl79l) which are 
adapted to the corners of our polytope P: 

Vp,k,s,{x)= n '^h'^s^^^A^)) (129) 

and 

d 

Vp,,{x) = \[k-\27,l,{x)fl'V^k^^{M,{x)l (130) 
i=i 

When we straighten out the corners by affine maps to be standard octants and separate 
variables x = {x',x") into directions near and far from dP, then Vp^k,Sk{x) is by definition a 
function of the near variables x' and Q^^s^ix) is by definition a function of the far variables 
x". 

The main result of this section is: 
Proposition 6.1. We have 



Phk{a) = Cmk"^ J det G^i^)Vp4^) (l + i?fc(p /i)) , (131) 
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where = 0{k~^) and Cm is a positive constant depending only on m. The expansion is 
uniform in the metric h and may be differentiated in the metric parameter h twice with a 
remainder of the same order. 

Equivalently, with defined in /il26\) and by letting 6k = k~s ^ 



VAc^) = C^A:^(--^*(^»y^^I(|)Pp,M.(f) {l + Rk{^,h)), (132) 
where again = 0{k~^). 



The factor k^ ("^~^k (f )) is due to the fact that we apply complex stationary phase in 
m — 6f{^) variables to a complex oscillatory integral with symbol of order fc(™-~^fc (f 

As a check, let us consider the m-dimensional Bargmann-Fock case where = r, and 

with no loss of generality we will assume that the first r facets are the close ones. The factor 
k"^ in the symbol of the Szego kernel is then split into k^ (absorbed in Vp^k,Sk] 
the far factor. As discussed in §3.0. ![ the far factor should have the form 




j=r+l j=r+l 

The asymptotic factor in Proposition 16. ![ 

n 

j=r+l 

matches this expression. Here, and throughout the proof, we always straighten out the corner 
to a standard octant when doing calculations in coordinates. 

Secondly, as a check on the remainder, we note that it arises from two sources. As will be 
seen in the proof, in 'far directions' the stationary phase remainder has the form 0{ ^^^} Qp^ ) 

while in the near directions it has the form 0{k(d{j, dP))"^). When d{^,dP) ~ k~^ the 
remainders match. 

We break up the proof into cases according to the distance of ^ to the various facets as 
k oo. Since we are studying joint asymptotics in (a, k), a may change with k. 



6.2. Interior asymptotics. 

• — is (5-far from all facets 

k 

We first consider the case where ^ is 5-far from all facets as an introduction to the problems 
we face. In this case, we obtain asymptotics of the integral flll4p by a complex stationary 
phase argument. But it is not quite standard even in this interior case. In the next section, 
we goo on to consider the same expansion when 6 depends on k. 

Lemma 6.2. Assume that there exists S > such that ^j(f) > S for all j , i.e., that ^ is 
5- far from all facets. Then there exist bounded smooth functions A_j(x) on P such that 

V,.{a) ~ Cmk"^ ^JdetG.i^) (l + + + ■■■ + Os{k~'')^ . 
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Here, G^p = V^m ( ^2.^) and G^{^) is its value at p- its norm is 0{S ^) and its determinant 
zs 0(5-'"). 

Before going into the proof, we note that the only assumption on the limit points of ^ 
is that they are 6-fai from facets. The lattice points a are implicitly allowed to vary with 
k. Asymptotics of the left side clearly depend on the asymptotics of the points ^, and the 
Lemma states how they do so. 

Proof. We now apply the complex stationary phase method, or more precisely its proof. The 
usual complex stationary phase theorem applies to exponents k^{6) where ^{0) is a positive 
phase function with a non-degenerate critical point at 6' = 0. In our case, the phase is also 
/c-dependent since it depends on ^ and the asymptotics of fll20p therefore depend on the 
asymptotics of ^ in the domain d{^,dP) > 6. Our stated asymptotics also depend on the 
behavior of ^ in the same way. 

Although the exact statement of complex stationary phase |Hoj (Theorem 7.7.5) does not 
apply, the proof applies without difficulty in this region. Namely, we introduce a cutoff 
Xs{(^) = x(^~^^) ^ (^ooj^x™) which = 1 in a 5-neighborhood of ^ = and which vanishes 
outside a 25- neighborhood of ^ = 0. We decompose the integral into its xs and 1 — Xs parts. 
A standard integration by parts argument, essentially the same as in Lemma [1.21 shows that 
the 1 — X(5 term is = 0{5~^^k^^) for all M > 0. In the x& part the integral may be viewed as 
an integral over and we may apply the Plancherel theorem as in the standard stationary 
phase argument to obtain 

(133) 

where Tg^^ is the Fourier transform. 



The stationary phase expansion (see [Ho] . Theorem 7.7.5) is the following: 



k^'v^,Mf^-,\T)^ ^^V.-^(7:), 0) |.=o (134) 



f \m/2 

'idetifs • 

k I 



where 



P^,,Afc(0)= V y^'-^{H-^^De,DgY{A^R^)\g=o (135) 

u—fi=j 2i'>3fi 

The only change in the standard argument is that we have a family of quadratic forms 
Ha depending on parameters {a, k) rather than a fixed one. But the standard proof is valid 
for this modification. As in the standard proof, we expand the exponential in fll33p and 
evaluate the terms and the remainder of the exponential factor just as in [Ho] Theorem 
7.7.5, to obtain fll34p . which becomes 

(det(r ^G,(f ))) E-io {{G.mDe, Do))' x,e'=^3(M,.-Hf (;.-i(f ), e^V.-^(f ), 0, A:) 
-FO(A;-^^sup,,,„,,,^ |(0^(f)(D,,D,)^ X.e'=«^('=''''^"(^»^.(/^;:'(f),e^V.-'(f),0, A:)) . 

(136) 

Here, Gx is the Hessian of the symplectic potential, i.e., the inverse of H^-ii^^y (cf. 12. 2p . We 
recall that G^ has poles of order one when x G dP. When d{^,dP) > 6, its norm is 



0=0 
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therefore 0{S~^) and its determinant is 0(5"™). Since R3 vanishes to order 3 at the critical 
point, the terms of the expansion can be arranged into terms of descending order as in the 
standard proof. If we recall that the leading term of S is k"^, we obtain the statement of 
Proposition 16.11 in the 5-interior case. 

□ 

a 

• — is 6k-isLr from facets with k6k 00 

k 

We continue to study the complex oscillatory integral fllMp but now allow ^ to become 
(5fc-close to some facet, and obtain a stationary phase expansion (with very possibly slow 
decrease in the steps) under the condition that k5k 00. This should be feasible since 
the phase /c$| is still rapidly oscillating in this region, albeit at different rates in different 
directions according to the proximity of ^ to a particular facet. The principal complication 
is as as follows: 

• The Hessian G^(^) now has components which blow up like 5^^ near the close facets. 
In the stationary phase expansion, we get factors of 

both in the expansion and remainder. We must verify that these terms still are of 
descending order. 

As a guide, we note that by fl^JS]) . the Bargmann-Fock phase with fi{z) = f is given by 
<^BF,f{0) = e'' - ^0) = {cose + z(sin^ - 6), ^), 

while the amplitude is constant. In this case, the phase factors into single-variable factors 
and one can employ the complex stationary phase method separately to each. In the general 
case, we will roughly split the variables 6 into two groups [6', 6"), depending on ^, so that 
the 6' variables are paired with the small components of ^ while the 6" variables are paired 
with its large components. The complex stationary phase method applies equally to either 
dd' or d6" integral, but the orders of the terms are determined by the proximity of ^ to the 
facets. 

Lemma 6.3. Let {5k} be a sequence such that kSk 00. Assume that > Sj. for all j, 
i.e.,, that ^ is 6k far from all facets. Then in the notation of Lemma \6.S\, we have 

VAo^) - Cr^k^ ydetG,(^) (1 + + + ■ ■ ■ + + 0(M,)-^^) , 

where now 

A^,{^)<D6-' = Cdq,dPr. 

Remark: One may regard this as an expansion in the semi-classical parameter {k6k)~^ = 
{kd{%dP))-\ 

Proof. We need to prove that the expansion fll57p may be re-arranged into terms of decreasing 
order and that the remainder can be made to have an arbitrarily small order k~^'^ by taking 
sufficiently many terms. 
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To analyze the expansion (11571) . we begin with a decomposition of the inverse Hessian 
which is the Hessian of the symplectic potential, which has the form uo + g where g G C°°{P) 
and where uq is the canonical symplectic potential fl56|) . We continue to fix a small 6 > 
as in the previous section, and consider the facets to which ^ is (5-close. We use the affine 
transformation to map these (5-close facets to the hyperplanes Xj = 0. In these coordinates, 
we may write the symplectic potential as 

"v(^) = X] 4 + gix), (137) 

where the Hessian of g is bounded with bounded derivatives near f . The Hessian Gs. then 
decomposes into the sum, 

GA^)= E \Sn+'^'9--=Gi + V'g, (138) 

where V'^g is smooth up to the boundary in a neighborhood of J^s^i^)- The notation G'^(x)^ 
refers to the 'singular part' of Gx- The choice of 6 is not important; we are allowing ^ to 
become 6k close to some facets, and for any choice of S, the sum will include such facets. 

The decomposition f ll38p of the inverse Hessian induces a block decomposition of the 
Hessian operator {GslDq^Dq). The chage of variables to x above induces an affine change 
of the 9 variables, as follows: We are using the coordinates [x' ,x") on P with x' denoting 
the linear coordinates in the directions of the normals to the facets The normals 

corresponding to ^^(f ) generate the isotropy algebra of the sub-torus (T™)' fixing the near 
facets. We have T™- = (T™)' x (T™)", and denote the corresponding coordinates by [6' , 9"). 

The Hessian operator in these coordinates has the form 

{G^i^)De,De)= Yl ^ D^,,,, + {G^{^)" De, De) , (139) 

where the second term has bounded coefficients. Evidently, the change to the interior sta- 
tionary phase expansion is entirely due to the singular part of the Hessian operator, 

{G'^De,De):= Yl ^ D^,,. (140) 

We now consider the order of magnitude of the terms in the jth term fll35p . which has 
the form 

k-" i{G^{^)De, DerAk{fi-,\^), e^V.-'(f ), 0, k)R,{k, 9, f^-'i^W) \e=o (Ml) 

with u — fi = j and with 2z/ > 3fi. The latter constraint is evident from the fact that R-s 
vanishes to order 3. 

Using (11391) . {Gip{j)D0, DqY becomes a sum of terms of which the most singular is 

(G;(f)D,,D,r:=( KDi,,r- 
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We will only discuss the terms generated by this operator; the discussion is similar but 
simpler for the other terms. In the extreme case of {G'LDg, Dg)'^, the discussion is essentially 

the same as in the previous section; in particular, (11411) has order . 

The problem with each application of {G^^{^)D0, Dg) is that it raises the order by the 
maximum of which may be as large as k6k- Although we have an overall and con- 

straints v — fi = j,2u > 3/i, it is not hard to check that these are not sufficient to produce 
negative exponents of k. 

The key fact which saves the situation is that the phase and amplitude S depend 
on 6 as functions of e*^|/x~^(^)p. Although has a more complicated 6'-dependence, its 
third and higher derivatives are the same as those of $|, and it is obvious that only these 
contibute to (11411) . Hence derivatives in 6 bring in factors of |Ai~^(f)P by the chain rule. 
Due to the behavior of the moment map near a facet, these chain rule factors cancel a square 
root of the blowing up factor in G^{^). This turns out to be sufficient for a descending series 
due to the power k~^ and constraint 2z/ > 3z/. 

Before giving all the details, let us consider what should be the 'worst' terms of (11411) . i.e., 
the ones with the least decay in k. Each factor of R3 comes with a factor of k, so one would 
expect terms with large fi to be 'worst'. The 'worst' term will be one with a maximum fi and 
where a maximum number of applications on operator {G%_Dg, DgY is applied to the 'chain- 
rule' factors (e*^|Ai"^(f ))jP (the jth component of this vector), obtained from an application 
of some Dgi to S or to R^. If instead we differentiate 5* or R^ again, we pull out another 

chain rule factor, which cancels more of the bad coefficient \. 

We now give the rigorous argument. The terms of (I14ip have the form, 

k-''^^' GJ'^Y'^' ■ ■ ■ GJ^y^^^D^'Rs ■ ■ ■ Df^^R^D^^+'S, (142) 
k k 

where = 2z/ and where D^'^ denote universal constant multiples of the multinomial differ- 
ential operators Qgr, f'^lg>'i3q where the union of the indices agrees with {^i, ji, . . . , iu,ju}- We 
need each \l3q\ > 3 for q < fi to remove the zero of -R3. If we only consider the most singular 
term, then we need iq = jq E In this case our term becomes 

{ n -^^D^'R-i-'-D^^^R^D^^+^S, (143) 
V=i;9.e.F.(f ) ""^ / 

For each factor there exist two factors of the associated differential operator When 

one is applied to either i?3 or S it pulls out a chain rule factor e^''^ |/i~^(^))jp. The second 
could be applied to this factor, hence need not introduce any new factors of |/^~^(f ))jP- We 
now estimate ( I143P by 

Now /i~^(a;) = Vu^{x) in p coordinates. So the square of the g^-th component of /i~^(|) 
equals log ^ plus a bounded remainder in p coordinates; here as above we are using the Xj 
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coordinates adapted to p It follows that in the z coordinates adapted to the facets of V 
corresponding to the hyperplanes x'^ = 0, with \zj\'^ = e^^, < C^. The constant 

C comes from the smooth part of the symplectic potential and has a uniform bound. As 
a check on the square root, we note that for the approximating Bargmann-Fock model we 
have \zj\^ = ^. It follows from I^SM and ^ < Cd{^,dP)-^ that 



«l < Ck-^^' (n;=i;,,e^,(|) 4) UU 4- ^ Ck-^^'d{f, dP)-'^' ^^^^^ 

= C{kd{%dP))-^, 

Effectively, the 'semi-classical parameter' has changed from k~^ to k~^d{^,dP)~^, a nat- 
ural parameter in boundary problems. As long as kd{^,dP) ^ 00 at some fixed rate, we 
obtain a descending expansion. 

□ 

6.3. Boundary zones: Corner zone . Having dealt with the case where |^| > 5^, we now 
turn to the complementary cases where d{fi{z), dP) < 6k, i.e., at least for one j, |^| < Sk or 
equivalently, ^ is (5fc-close to at least one facet. The choice of the scale 6k is so that it it is 
small enough to justify the Bargmann-Fock approximation in the 'near' variables. 

In this section, we consider the extreme 'corner' case where fi{z) lies in a (5fc-corner, i.e., 
where there exists a vertex v G dP so that d{fi{z),v) < 5k- Putting = 0, the assumption 
becomes that \^i{z)\ < C6k ■ Our main object is to determine the scale Sk so that the 
Bargmann-Fock approximation is valid. That is, for z = we should have in the 

multi-index notation of §2.61 (see flllSp ). 

~ P,|^(a) = A;-(27r)--/^„e^KE™i^^..^(^"^-i+^^.'^))d^. (146) 
Lemma 6.4. If fi{z) lies in a 6k-corner, then 

P,.{a) = C^P,|^(a)(l + 0(4) + Oik6l)) = C„P,.^(«)(1 + 0{k5l)). 
Proof. We may assume that v = and that the corner is a standard octant. The phase is 

k[Fc{\z\'e^')-F{\z\')-{^,e)). (147) 

We Taylor expand F{w) aX w = 0: 

Fc{e''\z\^) = F(0) + F'iO)e''\z\^ + Oi\z\^), 

so that 

Fc(|2|V^) - Fd^n = F'(0)|^|2(e^^ - 1)) + 0(1^1^). 

Since |2p = 0{6k), it follows that k times the quart ic remainder is 0{k6l) = o(l) as long 
as 6k = o(-^). Hence this part of the exponential is a symbol of order zero and may be 

absorbed into the amplitude. Further we note that F'(0)|2p = ^{z) + 0(|2;|^) and therefore 
we have 



k i^Fc{\z\'e'') - F{\z\') - t{-,e) ) = A;/i(z)((l - cos^) + z(sin^ - 6)) + 0(1^1"^)). 
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It follows that when fi{z) = f = 0{6k), the phase equals 

a{{l - cose) + i{sme - 6)) + 0{kSl). 

Absorbing the e^'^'^^^) = 1 + 0{k6l) term into the amplitude produces an oscillatory integral 
with the same phase function as for the Bargmann-Fock kernel. 

Now let us consider the amplitude of the integral. We continue to use the notation of 
Proposition 14. 2[ The amplitude has a semi-classical expansion Ak{z,w) ~ k"^ao{z,w) + 
k"^~^ai{z,w) + ■ ■ ■ . Further, the T™'-invariance implies that Ak{e'^^z,e^^w) = Ak{z,w). The 
leading order amplitude equals 1 when z = w and thus 

ao{z, e'^w) = 1 + Ce^^l^p + 0{\z\'^), 

hence the full symbol satisfies 

Ak{z, e''z) = + Ce''\z\^ + ■■■) + 0(61). 

When fj,{z) = f = 0{6k) we thus have 

A,{z,e^'z) = k"^ (l + Ce'^^ + 0{6l)). 

Therefore, = Vi,k («)(! + 0{6k) + 0{kSl)) in the corner region. 



□ 



6.4. Boundary zones: Mixed boundary zone . Now let us consider the general case 
where d{iJ,{z),dP) < 5k, but where is not necessarily in a corner. Thus, at least one 
component ^ = 0[5k) but not all components need to satisfy this condition. We refer to 
this CclSG clS 'mixed' since some components are small and some are not. 

The basic idea to handle this case is to split the components into 'near' and 'far' parts, 
to use Taylor expansions and Bargmann-Fock approximations in the near components, and 
to use complex stationary phase in the far components. By §6.2[ complex stationary phase 
works for any sequence 6k satisfying k6k — * oo, and by §6.31 the Taylor-Bargmann-Fock 
approximation works whenever 5^ = o(-^), so we have some flexibility in choosing 5k- 

Remark: In fact, we see that both the complex stationary phase and the Bargmann-Fock 
approximations are valid for k satisfying (for instance) ^ 5k < ^' ^logk ^ although the 

remainder estimates will not be equally sharp by both methods. In fact, the stationary phase 
remainder is of order {k5k)^^ while the Bargmann-Fock remainder is of order k5l] the two 

2 

remainders agree when 5k = k~3 and for small 5k the Bargmann-Fock remainder is smaller. 

We first choose linear coordinates so that fi{z) = | is 5^ close to the first r facets and 
5'f^ far from the p := m — r remaining facets, and by an affine map we position the first 
r facets as the first r coordinate hyperplanes at x = 0, and the remaining facets as the 
remaining coordinate hyperplanes. We use coordinates {x',x") relative to this splitting. We 
also write the z variables as {z', z") in the corresponding slice-orbit coordinates and (6'', 9") 
as the associated coordinates on T"^. 

We now introduce two small scales, a smaller one 5^ to define the nearest facets, and a 
larger one 5'^. The Bargmann-Fock approximation will be used in x' variables which are 5^ 
close to a facet. It is sometimes advantageous to use the Bargmann-Fock approximation also 
x" which are 5^' small, but the complex stationary phased method is also applicable. In the 
following, we continue to use the notation above Proposition 16.11 
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Lemma 6.5. Assume fi{z) lies in the mixed boundary zone {\x'\ < S'^, \x"\ < S'^}. If 
then Vhk^a) has an asymptotic expansion 

Our strategy for obtaining asymptotics of V^k (a) in this case is as follows: 

• We employ steepest descent in the p directions which are ^^-far from all facets, i.e., in 
the x" variables. This removes the x" variables and produces an expansion analogous 
to that of Lemma 16.21 

• In the remaining x' variables, we Taylor expand the phase and amplitude in the 
directions (5fc-close to dP as in §6.31 

• We thus obtain universal asymptotics to leading order depending only on the number 
of facets to which ^ is ^^-close. 

Proof. We are still working on the oscillatory integral with phase (11141) . but we now treat 
it as an iterated complex oscillatory integral in the variables {6', 6") defined above. We first 
consider the d6" integral, 

(148) 

where p is the number of 6" variables. We also let r = m — p be the number of 6' variables. 
We now verify that we may apply the complex stationary phase method to the d6" integral 
for fixed 9'. Throughout this section, we put z = yU~^(^) and often write Iz'l"^, {z"]"^ for the 
modulus square of the associated complex coordinate components of this point in the open 
orbit. 

The first step is to simplify the complex phase. As in §6.31 we Taylor expand -Fc(e*^ \z'\'^,e^^ 
in the z' variable (and only in the z' variable) to obtain 

Fc(e*^>f , e^^'V'f ) = Fc(0, e^^'V'f ) + F{(0, e*^'>'f )e-'^>f + 0(|z'|^), 

where Fi is the ^'-derivative of F. The phase is then 

k (Fc(e^^>f , e''"\z"\^) - F{\z'\^, \z"\^) - 6')) - z(f , 6")) 

= k (Fc(0,e*^'>'f ) - F(0, l^'f ) + k (Fi'(0, e''^"\z"\^)e'^'\z'\^ - F{(0, |^'f )|^f ) (149) 

-k {z{^,e') + z{^,e"))+o{k\z'\'). 

We now absorb the exponentials of the terms kO{\z'\'^), ki{^,6') of the phase (11491) into 
the amplitude, i.e., we take the new amplitude A'l to be the old one Ak multiplied by this 
factor. The term /cO(|z'|'^) is o(l), while ki{j^, 9') is constant in 9" , so their exponentials are 
symbols in 9" and may be absorbed into the amplitude. Morevoer, the term — F^(0, |2;"p)|2:'p 
is independent of 9" so its exponential may also be absorbed into the amplitude 

The phase function for the d9" integral thus simplifies to 

A;(Fc(0,e^^>'f)-F(0,|/f) +A;(Fi'(0,e^^'>'f)e*^>f) -A;z(^,0 (150) 
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Due to the presence of \z'\'^, the terms fc(F/(0, e'^" \z"\'^)e^^' \z'\'^ - F[{<d, \z2\'^)\z'\'^) are O(M'), 
hence of much lower order than the remaining terms. To simplify the phase further, we now 
argue that their exponentials can also be absorbed into the amplitude, albeit as exponentially 
growing rather than polynomially growing factors in k. Since -^{(O, |2;2p)|2'P) independent 
of 6'", it can be factored out of the 6" integral, so the key factor is 

E,(^"):=e'=(^»o,e-"|."P)e-'|.'P)^ (15^) 

where in the notation for Ek we omit its dependence on the parameters \z'\'^,9' . Thus we 
would like to show that complex stationary phase method applies to the complex oscillatory 
integral with phase 

:= (Fc(0,e'^'>'f ) - F(0, ) - (152) 
and with the amplitude A'l^{6") given by the original amplitude multiplied by 

expfc(^F{(0,e^^>'f )e^^>f -F{(0,|/f)|zf + z(^,0 + C'(A;|z'|^^ . 

The 'amplitude' is of exponential growth but its growth is of strictly lower exponential growth 
than the 'phase' factor. 

The next (not very important) observation is that by (11161) . the real part of complex phase 
damps the integral so that the integrand is negligible on the complement of the set 

modulo rapidly decaying errors. This follows by splitting up the sum in (ll22p -( fT23l) into the 
close facets to z and the far facets. The integrand is neglig ible unless \m\ < C^^; hence 
it is negligible unless 

d{e^'z,zf ~ E,e^,jM^))(l-^°«^")^;'(/iW)+0(kf) 

- E,.^.jM.))(^")^^;'(/^(^)) + 0(5D (154) 

Under the assumption that d" {jji{z),dP) > 5'1, the integrand is rapidly decaying unless 
^1 < ^F'' W^ could introduce a cutoff of the form x(-y/^^)) but for our purposes, it suffices 
to use a smooth cutoff x<5(^") around 9" = with a fixed small 6 so that we may use local 
6" coordinates. We then break up the integral using 1 = + (1 — xs)- The (1 — xs) term 
is rapidly decaying and may be neglected. 

We observe that V0"Fc{O,e^^ k"P) = "^/^cd-^"!' \^"\) is the complexified moment map 
for the subtoric variety z' = 0, and we can use Proposition 12.11 to see that its only critical 
point in the domain of integration is at 6" = 0. We denote the Hessian of the phase (11521) 
at 6" = by 

i/|':„|. = Vl,<^"{e")\e"=o = V^„Fc(0,e*^>'f )|,»=o, (155) 
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and observe that it equals iDij!^{\z"\,e^^" \z"\), the derivative of the moment map from the 
subtoric variety to its polytope. By the same calculation that led to fll38p . the 9" — 9" block 
of the inverse Hessian operator has the form 

= E ^ + V'^? := Gl„ + y'g, (156) 

3=1 

where \x"\ > 5^'. 

We now must verify that the complex stationary phase expansion 

{deik-'G'!^{\zr)Y'"Y.f=,k-^ {{G^^^^^^^^ (157) 
is a descending expansion in well-defined steps and that the remainder 

k-^' sup |(G';(|z'f)(Z^,„,D,.)''A'fc(^")"'x5^fc(^")|- (158) 

9" i^suppxs 

is of arbitrarily small order as M increases. 

We first note that the Hessian operator k~^ {G'^{\z"\'^){Doii , Dqh) brings in a net order of 
^~^(^fc)~^' since the coefficients ^ in the singular part are bounded by The maximal 

order terms arise from applying the Hessian operator to the factor E^. Each derivative can 
bring down a factor of A;Fi'(0,e^^'>"ne^^>'n = 0(MX). Since there are two 9" derivatives 
for each ^"^(5^)"^ maximum order in k from a single factor of k~^ {G'^{\z"\'^){Dqii ^ Dqu) 
applied to A'^ is of order 

Vk = k-\6l)-\{k5',n6lf + k6',6l) = k{6i)X + Si- 

In particular this is the order of magnitude of the sub-dominant term. Therefore, to obtain a 
descending expansion in steps of at least k~'"^, we obtain the following necessary and sufficient 
condition on {S[,S'^): 

Vk < Gk-'°. (159) 

Under this condition, the series and remainder will go down in steps of k'"^". 

With these choices of (5^, 6'^.), the complex stationary phase expansion gives an asymptotic 
expansion in powers of k~'^°. Recalling that the unique critical point occurs at 9" = 0, the 
remaining d9' integral is given by the dimensional constant Cm(27r)~'' times 

n.(«) ~ (det(A;-iG";(|^'f ))'/V^.e^^a'^')ESifc-^' ((G'^d^'f )(^.", ^."))' X5^fc(^', 0)^^', 

(160) 

plus the integral of the remainder (11581) . which is uniform in 9' and integrates to a remainder 
of the same order. Here we wrote the amplitude as A'l{9',9") and set 9" = after the 
differentiations. 

The differentiations leave the factor Ek (11511) while bringing down polynomials in the 
derivatives of its phase. The same is true of the factor e'^'^^"^'" •* that we absorbed into the 
amplitude. We now collect these factors and note that the exponent is simply the original 
phase (11491) evaluated at 9" = 0: 

^'{9'; \z'\^, \z"\^) := Fc(e*^>'|2, \z"\^) - F{\z'\^ \z"\^) - 9')) (161) 
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We also collect the derivatives of this phase and the other factors of Ak and find that 

M 

where Ak{6') is a classical symbol in k whose order is the order m of the original symbol A^. 
The integral (11601) then takes the form 

Vh^ia) ~ C„(det(fc-iG';(|^'f)))'^'/T^e'=*'^''''''''''"''^^^(^')rf^'- (163) 

This is a corner type integral as studied in §6.31 with \z"\'^ as an additional parameter. The 
asymptotics of ( 11630 are given by Lemma [6741 It is only necessary to keep track of the powers 
of l-z'P, \z"\'^ and of the parameter k~^{6'^)"^{k6'i^)'^ in the analysis of Ak- 
To do so, we first observe that 

Ve'Fc{e^''\zf, |z'f ) = t^i'M^',z"), {e^'' z' , z'% (164) 

i.e., it is the ' component of the complexified moment map. By definition of {z', z") it equals 
Y when 6' = 0. It follows that F{(0, |z"p)|2'p = ^, and the almost analytic extension 
satisfies 

F{(0,|z'f)e^^>f = e^^'^, (165) 

k 

where (as previously) the multiplication is componentwise. If we then Taylor expand the 
phase, we obtain 

$'(^'; Izf , b'f ) = F{(0, \z"\^)\z'\^(l - e''') + 0(\z'\^) = ^(1 - e''') + 0(\z'\^). (166) 

k 

If we absorb the e'^'^*^'^'*'' factor into the amplitude, he integral has now been converted to 
the form fll46p with a more complicated amplitude. 
We next observe that 

Ak = k"'{l + 0{\z'\')). (167) 
Hence, the assumption {z'l"^ = 0(5^) implies that to leading order 

P/,fe(a) ~ ^detk-^ G'^{\z''\^)k"^ f^,.e-'i^''''-'-''^)'^d9'{l + 0{6[)) 

(168) 

= A;™-fy/d^tG^[^)n.^(a')(l + 0(510). 
This completes the proof of the Lemma. 

□ 



6.5. Completion of proof of Proposition [6711 We now complete the proof of Proposition 
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6.5.1. Asymptotic expansion for Vhk{a). The error terms for the asymptotics of Vhk {a), 
in the corner zone, the interior zone and the mixed zone, are given by k~^{6'j^)~^, ^(^I-)^ 
Vk = k~^{6'l)~^ + + k{6'i^)'^6'l + 6[ respectively. In order to minimize these terms, we 

let 

By elementary calculation, the optimal choice for and is given by 

6', = 6'l = k-'^ and ^"^(5^"' = M^D' 

and 

2 

We let 6k = k~3 and break up the estimate into four cases. 

(1) \x'\, \x"\ < 6k'- this is the corner case handled in Lemma [6.41 if k{6kY ~^ 0. 

P,.(a) = C^n.^(a)(l+0(A;"^)). 

(2) \x"\ > 6k- By Lemma [6.3[ stationary phase is valid and 



P,.(a)~C„A;f ^detG^(^) (l + 0(A;-i 
(3) < 6k and \x"\ > 6k- By Lemma 1631 



Vh^{a) = C™A;-"Sy/detG";(f)Pp,fc,5^(a')(l +0(A;-3)). 

(4) \x"\ < 6k and \x'\ > 6k- This case is the same as case (3) by switching x' and x" - 
Combining the formulas above, the asymptotics for V^k (a) is given by fll32p 



Vhk{a) 



where Rk{%h) = 0(A;"3). 

On the other hand, equation (11311) is derived by the following calculation. 



ki.C^) 



Uf .(a-]kej{f) 



= kf^detG^{^)Vp,kil + 0{k-'^)), 
where the last equality follows from the Stirling approximation. 
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6.5.2. Derivatives with respect to metric parameters . Now suppose that h = ht is a. 
smooth one-parameter family of metrics. We would like to obtain asymptotics {■§^yVi^k(a) 
for j = 1,2. 

Proposition 6.6. For j = 1,2, there exist amplitudes Sj of order zero such that 

where Rk = 0{k~^). The expansion is uniform in h and may be differentiated in h twice 
with a remainder of the same order. 

Proof. Such time derivatives may also be represented in the form (11141) 

Ob J 'pm 

with a new amplitude A^j that is obtained by a combination of differentiations of the original 
amplitude in t and of multiplications by t derivatives of the phase. It is easy to see that 
t derivatives of the amplitude do not change the estimates above since they do not change 
the order in growth in k of the amplitude. However, t derivatives of the phase bring down 
factors k{-^y {Ft{e^^ — Ft(/i^^^(^)). The second derivative can bring down two factors 
with j = 1 or one factor with j = 2. We now verify that, despite the extra factor of k, the 
new oscillatory integral still satisfies the same estimates as before. 

The key point is that, by the calculation flllSp . the phase Ft(e*^/i^^^(^)) — Ft{ixJ^^{^)) — 
i{^, 9) for any metric h vanishes to order two at the critical point 6* = 0; the first derivative 
vanishes because V eFi^e^^ z)\g=Q = ifih{z). Hence, the t derivative of the /it-dependent Taylor 
expansion f lllSp for a one-parameter family ht of metrics also vanishes to order 2, i.e.. 

Thus, for each new power of k one obtains by differentiating the phase factor in t one obtains 
a factor which vanishes to order two at 6* = 0. As a check, we note that in the Bargmann-Fock 
model, the phase has the form ^j(e*^^ — 1 — iOj)^. 

Let us first consider the first derivative. We repeat the asymptotic analysis but with the 
new amplitude 5*1. In the 'interior region' the stationary phase calculation in Proposition [^2] 
proceeds as before, but the leading term (now of one higher order than before) vanishes since 
it contains the value of fll7UI) at the critical point as a factor. Therefore the asymptotics start 
at the same order as before but with the value of the second ^-derivative of the amplitude 
at ^ = 0. 

In the corner, resp. mixed boundary, zone we obtain an integral of the same type as the 
ones studied in Lemma 16.41 resp. Lemma 16.51 but again with an amplitude of one higher 
order given by the t-derivative of the phase. The only change in the calculation is in the 
Taylor expansion of the amplitude in fll67p in the z' variable, which now has the form 
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SO that the final integral now has the form 

(2.)-"*-/^^e-'(<-''-'-"'))* (F,(e-V,;'(^)) - F'(P'.'(j))))^„_/0'- 
As noted in (117(1 

Ml) (i^*(e^Vr^(f )) - i^*(/ir^(f ))) = Ml) {Ft{e^'f^7\W - Ftif^T'm - ^(f ' ^)) 

= ^l/o(i--)£(^t(e^^Vr^(f))^^ 

= 0(A;|^pf). 

Since the stationary phase method applies as long as |a| — > oo we may assume that |a| < C 
and we see that the factor is then bounded. Here, we have suppressed the subscript C for 
the almost-analytic extension to simplify the writing. 

As an independent check, we use integration by parts in 6'. We use a cutoff function x 
supported near ^' = to decompose the integral into a term supported near 6' = and one 
supported away from 6' = 0. We use the integration by parts operator 

£ = Je''' -l)a'- Ve 

where we note that the factors of k cancel. The operator is well defined for 6' ^ and 
repeated partial integration gives decay in a' in case \a'\ — » oo. On the support of x the 
denominator is not well defined but the vanishing of the phase to order two shows that £*(S'i) 
is bounded. 

Now we consider second time derivatives. The second ^ could be applied to the phase 
factor e'^*' again or it could be applied again to fll7ip . and then we have 

(172) 

+ (Mft)(i^*(e^V.-^(f )) - Ftif.-,\fW + Oilzf). 

The first term contains the factor and after cancellation it induces a term of order ja'p. 
In addition this term vanishes to order four at 6' = 0. Hence the stationary phase calculation 
in the case of the first derivative equally shows that the first two terms vanish and thus the 
factors of fc^ are cancelled. In the regime where stationary phase is not applicable, may 
be assumed bounded, and additionally one can integrate by parts twice. Thus again this 
term is bounded. 

□ 

6.6. Completion of the proof of Lemma 11.31 So far we have only considered the asymp- 
totics of Vhk{t, z). We now take the ratios to complete the proof of Lemma [1.3[ 



Lemma 6.7. With 6^ defined by / f55]) . we have 

ltk[t,a) - (^(detV^no(f))i-'(detV^«i(f))*; [i- + ^ 

The asymptotic in may be differentiated twice with the same order of remainder. 
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Proof. Combining Corollary 13.21 and Proposition 16 . 11 we have 

n,{t, a) = v/detG,,(f) Pp,.(f) _ . 

(v/detG,„(f) Pp,.(f))i-*(v/detG,,(f) VpAW' 
We observe that the factors of Vp^k cancel out, leaving 

(vdetG^^d)! *(^detG^,(f))* 

By Proposition 16.61 the asymptotic in fll73p may be differentiated twice with the same 
order of remainder, completing the proof. 

□ 

Remark: By ( !58|) . we also have 



Indeed, the factors of are independent of the metrics and cancel out. Also 
is smooth on P. 

The following simpler estimate on logarithmic derivatives is sufficient for much of the proof 
of the main results: 

Lemma 6.8. We have: 

(1) dtlog7lk{t,a) is uniformly bounded. 

(2) df\og7lk(t,a) is uniformly bounded. 

Proof. We first note that 

9^1og7^fc(t,a) = logP^fe(a) - logP^fe(a) - 9* logP^fe(a). (175) 
We note that by Proposition l6.lt 

1 (1/ ~ (1/ 1 

logV^Aa) = -\ogdet{k-'GA-)) +logPp,fc(-) + logC^ + 0(A;-3). (176) 



As in Lemma 16.71 the Bargmann-Fock terms cancel between the Hq and hi terms, while 
the metric factors simplify asymptotically to | log {S^-^6^g), and this is clearly bounded. To 
complete the proof of (1), we need that the final ratio is bounded. By Lemma [6.61 we see 
that in the 'interior' region both numerator and denominator have asymptotics which differ 
only in the value of a zeroth order amplitude at ^ = and that it equals 1 in the case 
of the denominator. Hence, the ratio is bounded in the interior. Towards the boundary, 
the denominator is comparable with the Bargmann-Fock model and is bounded below by 
one. The numerator is also bounded by Lemma 16. 6[ and therefore the ratio is everywhere 
bounded. 

Now we consider (2), which simplifies to 

a,Mog7^,(^,a) = -4-4V+KT^ • (i77) 
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As we have just argued, the second factor is bounded. The same argument apphes to the 
first term by Lemma 16. 6[ 

□ 

7. AND C^-CONVERGENCE 

We begin with the rather simple proof of C°-convergence with remainder bounds. 
7.1. C°-convergence. 

Proposition 7.1. ilogZfc(t,^)e-'=^'(^) = O(^) uniformly for {t,z) e [0,1] x M 

The Proposition follows from the following: 
Lemma 7.2. (Upper/Lower hound Lemma) There exist C, c > so that 

c < nk{t,a) < C. 

Proof. This follows immediately from Lemma 16. 7[ □ 

C°-convergence is an immediate consequence of the upper and lower bound lemma: 
Proof. By the upper/lower bound lemma, there exist positive constants c, C > so that 

cli^.{z,z)< n,{t,a) ' <CU,.{z,z). (178) 

Hence, 

\Sc{z)\l^ 

\\ogYiy^u{z,z) < l\ogJ2aekPnz^'^k{t,a) ^ ' < i log 0^,^(2, 2) + 0(i) 

(179) 

= 0(if^). 



where the last estimate follows from (1^5]) . 



□ 



7.2. C^-convergence. We now discuss first derivatives in {t, z). In the z variable the vector 
fields ^ vanish on T), so can only use them to estimate norms in directions 5k far from the 
boundary. In directions close to the boundary we may choose coordinates so that derivatives 
in z' near z' = define the norm. 

The estimates in the p and z' derivatives are similar. We carry out the calculations in 
detail in the p variables and then indicate how to carry out the analogous estimates in the 
z variable. 

We also consider t derivative. The key distinction between t and z derivatives is the 
following: 

• z or p derivatives bring down derivatives of the phase, which have the form k{iit{z) — 
^). The factor of k raises the order of asymptotics while the factor {pt{z) — f ) lowers 
it by the Localization Lemma. 

• t derivatives do not apply to the phase and only differentiate 7^fc(t,a) and Q/ife(a). 

Proposition 7.3. Uniformly for {t, z) e [0, 1] x M, we have: 
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(1) I. 



dp, 



0{t 



(2) The same estimate is valid if we differentiate in jf- in directions nearV as in Propo- 



sition 14- 6[ 



(3) I 



GfcPnZ' 



TZkit,a 



\So.iz)\l, 



0{k-^. 



Proof. We first prove (1). 



Vplog ^ 7^fc(^ 
y 



a] 



Qy,k (a) 



E 



'}Zk{t,a] 



+s 



E 



t 

t 



+ Oik 



0(A; 



proving (1). In this estimate, we use the Locahzation Lemma [L2] and the upper/lower bound 
Lemma [7.21 on TZk. 

Regarding derivatives in (2), the only change to the argument is in summing only a 
with a„ 7^ and then changing a ^ a — (0, . . . , 1„, . . . , 0) as explained in Proposition 14.61 
Clearly the localization and the estimates only change by ^. 

We now consider the dt derivative. By Proposition 14.41 we have 
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1 d 



1^* 



E. 



E„7^fe(i,a)9tlog7efe(t,a)^ 



< a.P> 



E. log Q, (a) (t,a)9t log Q , (a) ^i^^ N 



\ 



E. 



e<".P> 



log 



;i80) 

Notice that Q,-^k = Tlk{t, a){Q,^k{a)y and so dt\.ogQk{t,a) = dt\ogTlk{t, a) + 

. It follows easily from the fact proved in Lemma 11.31 (or more precisely the 



simpler Lemma [678]) that TZk{t,a) = 0(1) and dtlog(Tlk(t,a)) = 0(1). Also log ^ l(^ta) 
0(A;) uniformly in a. Replacing TZi^ by TZ^o plus an error of order k~3, 



17)7 2^ 7^fc(^, 



|5'q(2) 



aekPni 



Qhk{a 



□ 



8. O^-CONVERGENCE 

We now consider second derivatives in p, t. Again we must separately consider derivatives 
in the interior and near the boundary. The following Proposition completes the proof of 
Theorem 11.11 

Proposition 8.1. Uniformly for (t, z) e [0, 1] x M, we have, for any 6 > 0, 



(2) I 

(3) I 



'}Zk{t,a 



dtdpj 



\sJ{z)\l^ 



0{k- 
0{t 



0{t 



(4) The same estimates are valid if we replace in directions near D as in Proposition 

We break up the proof into the four cases. To simplify the exposition, we introduce some 
new notation for localizing sums over lattice points. By the Localization Lemma 11.21 sums 
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over lattice points can be localized to a ball of radius 0{k^^^^) around jJt{z)- We emphasize 
that although there are three metrics at play, it is the metric ht along the Monge- Ampere 
geodesic that is used to localize the sum. We introduce a notation for localized sums over 
pairs of lattice points: let 

\f-f,t{z)\,\&~^t{z)\<k-i+' 

8.1. Second space derivatives in the interior. In this section we prove case (1). We 
have, 



dpidpj 



-log 7^fc(^, 



a] 



\Sa{z)\l. 



aekPnZ" 



(182) 



OpiOpj 



-T 

2 



[a 



=<Q,p> 

OT) 



modulo O(^) by Proposition 14.31 We also completed the square and used that the sum over 
a is a probability measure to replace — a[3 by |(a — ■ We also use Lemma fOl to write 

dp^dp '^^ ^ lattice points. 

By the Localization Lemma [1. 2 [ each sum over lattice points can be localized to a ball of 
radius 0(A;~2+^) around pt{.z). Then, by Lemma [T73l each occurrence of 7^fc(t,a) or 7lk(t,P) 
may be replaced by TZooit, |) plus an error of order k~3. Since ^(a — P)"^ = Oik"^^) the total 
error is of order /c^^~i. Since 5 is arbitrarily small, this term is decaying. Further, after 



replacing TZk{t,(3) by TZooit, ^) we may then replace ^, | by pt{z) at the expense of another 

ICC 



error of order k 2+*. By modifying (11821) accordingly, we have 



dpidp. 



■log ^ 'R'k{t,a) 



Q,fc(a) 



+ 0(A;-3+25) 



|E.,,(«-/5)^ ^°o(t,/i,(e-/2))2^^ |E.,,(«-/5) 



E.^oo(t,/i*(e''/2))- 



e(a.p> 



(183) 



= 0, 



where = means that the lines agree modulo errors of order 0{k~^~^'^^). In the last estimate, 
we use that TZooit, Pt{G^^'^)Y cancels out in the first term. This completes the proof in the 
spatial interior case. 
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The modifications when z is close to dP are just as in the case of the first derivatives. 
8.2. Mixed space-time derivatives. The mixed space-time derivative is given by 



92 



log ^'^•(^' 



a] 



^Ea,/3(«-/5) 7^fc(^,/3)7^fc(t,«)9^log 



dpidt 



It suffices to prove that 



Ea/3(«-/5) 5^log(7^fc(^,«)) 



and 



1 

k 



/3) 7^fc(^,/?)7^fc(^,a)9^log 



0(A;- 



The first estimate follows by the Localization Lemma [1.21 and from Lemma [6. 8 [ i.e., that 
9t log (7^fc(t, a)) = 0(1). The second estimate is very similar to that in §8.H specifically in 
(11831) . so we do not write it out in full. In outline, we first apply the Localization Lemma and 
replace each TZk{t^ a) by lZoo{pt{z)) with z = e^/2_ The errors in making these replacements 

are of order because dtlog (^-^-pa)^ = 0{k\ut{a) — Ut{P)\) = 0{k^~^^) and because 

{a — /3) = 0{k^'^^) in the localized sum. We then express -^^^t in terms of the Szego 
kernel, i.e., as a sum over lattice points, using Proposition 14.51 and cancel the '^t term. 
The sum of the remainders is then of order completing the proof in this mixed case. 

8.3. Second time derivatives. The proof in this case follows the same pattern, although 
the estimates are somewhat more involved. The main steps are to localize the sums over 
lattice points, to replace each TZk by TZoa-, then to cancel out T^oo after all replacements, and 
to see that the resulting lattice point sum cancels k ^jj^t- The complications are only due 
to the number of estimates that are required to justify the replacements. 
The second time derivative equals 
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1 92 



log ZlaGfcPnZ™ T^kit, <y) Q (^a) 



(184) 



v 



Here, we have simplified the numerator of the first term by replacing 



9f log 




which is valid since the expression is ant i- symmetric in and since we are summing in 

To simplify the notation, we now abbreviate TZ{a) = TZk{t,a), T[a) = q 1{q) j f ~ %. 

h. 



at ' 



TV 



and write fll84p - ^ 
TV 



where the numerator has the schematic form 



;i85) 

from 



and where the denominator is D = iJ2a'^i'^)'^i'^))^ ■ omit the factors ^ ^^^^ ^ ^^^-^ 

the notation since they are always present. 

We now compare and D to the corresponding expressions in the second time derivative 
of the Szego kernel in Proposition 14.51 In the latter case, 7^ = 1 so any terms with t- 
derivatives of TZ above do not occur in the third comparison expression of Proposition 14.51 
Terms with no t derivatives of TZ will be precisely as in the comparison except that TZ is 
replaced by 1. So we consider the sub-sum of iV, 

, 2^ 



(186) 



If we now replace all occurrences of TZk{t,a) by TZoo{fJ't{z)) in both numerator and denom- 
inator we get the Szego kernel expression (the third comparison expression of Proposition 
14. 5p of order (This is verified in more detail at the end of the proof). So we are left with 
estimating two remainder terms: First, the difference A^i — Ni where Ni is a sum of terms 
in which we replace at least one TZ{a) by TZooif^t{z)) (or with /5). Second, we must estimate 
N-Ni. 

We first consider Ni—Ni. It arises by substituting at least one 7l{a) —TZoo{fit{z)) = 0{k~s) 
for one of the 7^(a)'s in A'^i. We apply the localization argument Lemma [1.21 to replace Ni 



(and D) by sums over p| G B{nt{z),k 2+"^). We thus need to estimate the following 
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expression, when at least one 7l{a) is replaced by 71(a) — 7loo{fk{z)): 



E„«7^fc(^,/3)7^fc(^,a) L9jiog 



k 



= {/3) 



, k (a) 



-g,logQ,.(«)) 



E 



a{a,p) 



( 



1 



o{a,p) p{fi,p) 



\ 



\ 



Due to the factor ^ outside the sum, it suffices to prove that 



T' 



la 



1 /T 



2 V T 



0{k 



l+25\ 



By Proposition 13. ![ we have 
Since Ut = {1 — t)uo + tui, we have 



T 



(«) = + M«i - «o)(^) = + M/i - /o)(^), 



where we recall from §2.21 that = uq + with /(^ smooth up to the boundary of P. 
It follows that, 



V 



a. 



(«) + + M/i - /o)(t) - M/i - /o)(t) 



V 



k' 



k' 



with 



(Y(«)y = -(^y = o(i). 



M/i - /o)(^) - M/i - /o)(f ) = ^0(1^ - ^1) 



(187) 
(188) 



Further, by Lemma [6.71 (using Lemma [6.6p . the factors of 

(|)^/.?(«) _ (5i(t,a,A;) + i?fc(f,/i)) 



So{t,a,k) 



0{k 



0(1), 



and similarly (^)' = 0(1). Since fll87p is squared, it has terms as large as 0{k^'^ ). Taking 
into account the overall factor of \ and the presence of at least one factor of size k~^ coming 
from the replacement of at least one IZkii, a) by 7loo{fit{z)) , we see that Ni — Ni has order 
^-3+2(5 gj^^ again this decays for sufficiently small 6. 
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Now we estimate N — Ni, which consists of terms with at least one t-derivative of TZ. By 
Lemma I6.7[ the terms with no t derivatives on T give the terms 

2 zi f 

k / ^2 



(S,fc(a)) 



1 J:a,f^T^k{t, P)n,it, a)d^ log (7^,(t, a)) f^' f^' 
+ ^ '-^ = 0{k"'] 



(S^fc(")) 



by Lemma IL3[ 

This leaves us with the terms 

TZ' 7^' T T 

Again by Lemma |6.8[ the first term is 0(1) while the second factor is fll87p and has size 
kk~2~^^. Here, we again use Propositions 13.11 and 16.61 Due to the overall factor of ^ this 
term has size k~^~^^. 

Therefore, as stated above, up to errors of order k~^^^~^^, fll84p is simplified to — 
plus 



0(01, p) e(/3,p> 



(a.p) 



(189) 

As before, we cancel the factors of 7^oo(/^t(e''''^)). The resulting difference then cancels to 
order /c^^/^^'^ by Lemma [4.51 (3). 

This completes the proof of the second time derivative estimate, and hence of the main 
theorem. 
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